Accepted by ApJ 3/22/11 

Preprint typeset using I^^'T^]X style emulatcapj v. 8/13/10 



O 



(N 

o 

u 

o 



> 
m 

00 

rn 
O 



HYDRODYNAMIC MODELS OF RADIO GALAXY MORPHOLOGY: WINGED AND X-SHAPED SOURCES 

Edmund J. Hodges-Kluck^ & Christopher S. Reynolds^ 

Department of Astronomy, University of Maryland, College Park, MD 20742-2421, USA 

Accepted by ApJ 3/22/11 

ABSTRACT 

We present three-dimensional hydrodynamic models of radio galaxies interacting with initially re- 
laxed hot atmospheres and explore the significant off-axis radio lobe structures which result under 
certain conditions. With a focus on the "winged" and "X-shaped" radio galaxy population, we con- 
firm the importance of observed trends such as the connection of wing formation with jets co-aligned 
with the major axis of the surrounding atmosphere. These wings are formed substantially by the 
deflection of lobe plasma flowing back from the hot spots (backflow) and develop in two stages: su- 
personic expansion of an overpressured cocoon at early times followed by buoyant expansion at later 
times. We explore a limited parameter space of jet and atmosphere properties and find that the most 
prominent wings are produced when a decaying jet is injected into a small, dense, highly elliptical 
atmosphere. On the basis of this search, we argue that the deflection of backflow by gradients in the 
hot atmosphere is a strong candidate for forming observed wings but must work in tandem with some 
other mechanism for forming the initial wing channels. Our models indicate that lobe interaction with 
the hot atmosphere may play a dominant role in shaping the morphology of radio galaxies. 
Subject headings: galaxies: jets - galaxies: active - hydrodynamics - galaxies: clusters: intracluster 
medium 



1. INTRODUCTION 

Recently, X-shaped radio galaxies (XRGs) — a peculiar 
class of radio galaxies with two pairs of m isaligned lobes 
(jEkers et al.lll978l:lLeahv fc Parmalll992D — have received 
significant attention as new observations have explored 
and critiqued the competing models for the origin of their 
odd morphology. The distinctive "X" shape occurs due 
to the intersection of two centro-symmetric synchrotron- 
emitting structures at a common nucleus (e.g. Figured]). 
One of these structures is associated with an active rel- 
ativistic jet (the "primary" lobes) whereas the other is 
fainter, more ragged, and does not appear to harbor a 
jet (the "secondary" lobes or "wings"). The wings can 
be long, collimated, and symmetric about the nucleus, 
and may even exhibit Z-sha ped morphology of their own 
(|Gopal-Krishna et al.ll200l . The origin of the wings is 
not presently understood, but XRGs share other charac- 
teristics which indicate a common origin. 

X-shaped sources mak e up about 5-10% of Fa naroff- 
Riley Type I I CFR, II: iFanarofr fc RUctI [1971 radio 
galaxies (e.g. iLeahv fc Parmal 119921 ). The consensus 
sample is fairly small, bu t with the id entification of 
- 100 "candidate" XRGs (iCheund [20071) . a search for 
trends is now possible. Based on the combined con- 
sensus and candidate sample, XRGs tend to have radio 
powers and host optical luminosities near the FR I/II 
"brea k" in a Ledlow-Owen radio power- visual magni- 
tude (iLedlow fc Owed 19960 type plot (jCheung et al.l 
[200I ILandt et aLllmaT and reside in elliptical galax- 
ies with larger-than-average black hole masses (inferred 
from the M — a rel ation and broa.d lin es from the ac- 
tive galactic nucleus: [Mezcua et al.ll2010l ). XRGs do not 
seem to be in galaxies currently u ndergoing a merger 
based on inferred starburst history ( iMezcua "erall[2010l) 
and a spectroscopic search for broad emission lines and 
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dusty nuclei (jLandt et al.ll20ldl ). but the dynamic ages 
of their activ e lobes are you nger than the age of the 
starburst (Me zcua et"ani2010t ). The wings are preferen- 
tially CO- aligned with the minor axis of the host, whereas 
the prim ary lobes tend to be co-aligned with the ma- 
jor a xis ([Capetti et al.l [20021 : iSaripalli fc SubrahmanvanI 
120091) : this alignment has been confirmed in the X-rays 
which trace the distribution of hot gas in and around 
the host galaxies that provides the working surface which 
shapes the radio galaxy (Hodges-Kluck et al. 2010). The 
radiative age of the wings docs not seem to follow a clear 
trend, as some wings have steeper spectral indices than 
the primary lobes (im plying ageing) whereas others do 
not ()Lal fc Raoll2005() . Notably, some wings are actually 
longer than the primary lobes, and wings are only found 
in radio galaxies with strong bridges. 

A critical review of the observational data a nd ex - 
isting models is found in iGopal-Krishna et al.l ()2010D : 
here we briefly summarize the main threads. The ori- 
gin of the secondary lobes has been attributed to (i) 
a rapid reorientation of the spin axis of the supermas- 
sive black hole (SMBH) powering the jet (i.e. the 
wings are fossils), either due to a SMBH merger (e.g. 
Rottmann|[2001l:lZier fc Biermannl|200lHMerritt fc Ekers 
2002 ) or rapid precession (e.g. iDennett-Thorpe et al 
2002|) : (ii) redirection and collimation of "backflow" 



(spent jet plasma flowing back from the te r minal 
shocks) (H eahv fc WiUiams 1984; Worrall et all 119951 : 
iCapetti etal..,2002:, Kraft et al..2005^ : (iii) a binary AGN 
(|Lal fc Raoll2007D : and (iv) interaction of the jet with dis- 
turbed morphology (e.g. stellar s hells, p hase- wr apped 
remn ants of a merged disk galaxy; iGopal -Krishna et al.l 
[201^. Since XRGs are usually strongly bridged sources 
and are apparently aware of their environments, it is 
worth examining closely the hypothesis that the X- 
shaped morphology originates from an interaction be- 
tween the radio galaxy and its environment. In this pa- 
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Figure 1. 3C 403, a typical XRG showing long, faint wings 
emanating from the same core as the edge-brightened lobes 
(east -west axis). The red and green ellipses have the ellip- 
ticity and positio n angle of the host galaxy and ISM (from 
IKraft et al.ll2005l ') respectively (not to scale). 

per, we seek to test the viability of the backflow model 
with three-dimensional hydrodynamic simulations in el- 
liptical atmospheres. 



In the backflow model, wings are produced by a 
single AGN outburst with powerful jets as the back- 
flow is diverted along the steepest pressure gradient 
of the surrounding atmosphere (i.e. the minor axis). 
The wings either rise buoyantly or are driven in this 
direction (Lcahv & Wflhams 1984; Wor rall et all 119951 : 
IKraft et all 120051) or form as supersonic outflows along 
the direction of l east resistance from a ruptured over - 
pressured cocoon (jCapetti et al.ll200llZa^ni et al.ll2003l) . 
These scenarios naturally explain the observed correla- 
tion between the wings and the minor axis as well as the 
dearth of strictly FR I XRGs (the weaker FR I sources 
do not produce powerful back-flows). However, the back- 
flow model faces significant challenges, most notably the 
very long wings in many XRGs (often longer than the 
primary lobes). Whereas powerful jets are expected to 
drive through the surrounding medium supersonically, 
wings which expand buoyantly would do so at most tra n- 
sonically for most of their lives (Lea hv fc Parmalll992[ ). 

Because, in this model, fluid effects are primarily re- 
sponsible for the "X" shape, hydrodynamic simulations 
are an ideal proving ground. Although recent hydro- 
dynamic simulations of r adio galaxy morphology ex- 
ist in the literature, (e .g. Sutherland fc Bickncll' feoOTt 
iFalceta-Goncalves et al., 2010: Gaiblcr ct al. 2010), most 
are not concerned with the formation of lateral asymme- 
tries such as wings. In view of advances in simulating 
radio galaxies in the past decade, the mounting evidence 
that XRGs constitute a population demands a critical 
look at the backflow model with new simulations. In 



this paper, we discuss 3D hydrodynamic models of radio 
galaxies ignited in initially relaxed, elliptical atmospheres 
and examine how wings form. 

We begin by discussing our simulation setup and strat- 
egy for exploring wing production in Section 2, then 
present our model runs and describe the evolution of a 
generic winged source in Section 3. In particular, we find 
that favorable pressure gradients are necessary but not 
suflicient to produce an X-shaped source; the character 
and time evolution of the jet is equally important. In Sec- 
tion 4, we briefly discuss the missing and artificial physics 
in our simulations, and in Section 5 we assess the back- 
flow model in view of our results as well as discuss the 
implications for the broader picture of radio lobe mor- 
phology. Finally, in Section 6 we summarize our main 
results and conclusions. 

Throughout this paper we use the terms "secondary 
lobes" and "wings" interchangeably. We also use the 
term "winged" source to refer to any radio source with 
substantial symmetric off-axis distortions whereas "X- 
shaped" sources are a subset of winged sources with 
an axial wing-to-lobe length ratio of more than 0.8. 
This nomenclature reffects the view of the hydrodynamic 
model in which XRGs are indeed a subset of a broader 
category of distortions. Finally, we use the term "back- 
ffows" to refer specifically to actual fiuid fiows heading 
back to the nucleus from the jet heads, whereas we use 
"backfiow" to refer to the material in these fiows. 

2. HYDRODYNAMIC SIMULATIONS 

We use a parallelized version of the ZEUS code 
(jStone fc Normanlll992allbl ) for our hydrodynamic sim- 
ulations. ZEUS is a second-order (spatial accuracy) Eu- 
lerian finite-differencing magnetohydrodynamic (MHD) 
code which solves the standard equations of hydrody- 
namics. 



| + V.(pv)^0 

dv VP 

— + V Vv= V$ 

ot p 



V • (pve) = -FV • V, 



(1) 
(2) 

(3) 



for an ideal compressible fiuid and introduces artificial 
viscosity for shocks. The version we use, ZEUS-MPv2, 
is based on the National Center for Super commuting Ap- 
plications (NCSA) code described in Hav es et al.l (j2006). 
We use spherical polar coordinates (r,0,0) in the purely 
hydrodynamical mode for all our runs; the jets are there- 
fore injected from an inner boundary sphere with some 
small, but finite, radius rjnncr- We outline our basic sim- 
ulation setup below, followed by our strategy for explor- 
ing winged sources and a description of the evolution of 
a standard double-lobed source for comparison. 

2.1. Simulation Setup 

Hydrodynamic and MHD simulations of radio galax- 
ies interacting with their surroundings are common due 
to the possibility that energy deposited by the lobes is 
a mode of heating in AGN feedback scenarios. We do 
not incorporate relativistic jet physics since we are pri- 
marily concerned with lobe mixing and evolution, but 
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the robustness of this assuraption is explored in Sec- 
tion 4.1. We note here that non-relativistic, light, hy- 
personic hydrodynamic jets reproduce several essential 
features of jets. These include recoUimation shocks, ter- 
minal shocks at the jet head (associated with radio "hot 
spots" ) and back-flows of spent material from the termi- 
nal shocks which sheath the jet and produce lobes. Since 
we are interested in the evolution of back-flowing fluid 
in anisotropic environments, we adopt a non-relativistic 
purely hydrodynamic scheme. 

We inject the jets as bi-directional flows into an ini- 
tially hydrostatic, ellipsoidal atmosphere. The atmo- 
spheres we set up are ellipsoids with a polytropic equa- 
tion of state (7 = 5/3) and a 3D ;9 -model density profile 
(jCavaliere fc Fusco-Femian"olll976D given by 



potential $ is 



(1 + (r/ro)2)3/3/2 



(4) 



where /3 = 0.5, vq is the core radius, and po is the core 
density. To adapt this model for elliptical atmospheres, 
we make the radial density profile dependent on angle by 
adding a major-axis position angle PA and ellipticity e to 
the atmosphere, which we denote as for an elliptical 
atmosphere. To easily allow for triaxial atmospheres, 
we define the ellipsoid in Cartesian space (xe, Ue, Ze) 
and break up e along each axis, then transform back to 
spherical coordinates to obtain the radius using the 
grid coordinates (r, 6*, </>): 

_ ^xlil - e^f + yg(l - ey)^ + 

Tp — r {0) 



- r sin(0) 



ye = rsin(6') 



z^ = r cos{9) 



(1 ^max) 

cos(0)cos(PA) -I- sin(0) sin(PA) 



sin((/))cos(PA) — cos(0) sin(PA) 



cos(PA) -I- sin(PA) 



(6) 



(7) 



(8) 



where PA is nominally measured clockwise from the z- 
axis (in practice, it is only the difference between the jet 
axis and the major axis of the ellipsoid that matters). To 
normalize the size of the atmosphere, emax in Equation 5 
is defined as the largest value of e along each axis. For ex- 
ample, an atmosphere with large is elongated along the 
z-axis. Equations 4 and 5 may be understood as saying 
that introducing e changes the effective core radius along 
a given axis: ro,cff = fo(l ^ e„iax)/(l - e)- For example, 
if Cmax = i-z — 0.75, — 0.0, and tq = 1.0, then along 
the z-axis ro.oS = 1.0 and along the x-axis ro.oft — 0.25. 
This phrasing is convenient when comparing the major 
and minor axes rather than the quantities along the full 
range of angles and and we use it hereafter. 

The atmospheres are initially set up in hydrostatic 
equilibrium. 



V$ = -Vp, 
P 



(9) 



assuming that the background dark matter potential is 
static and dominant so that the gas self-gravity is not 
important. We define the adiabatic sound speed Cg = 1.0 
in code units and set up an isothermal atmosphere so the 



$ = ^ ln(p) 

7 



(10) 



Since XRGs are powerful radio galaxies and exist at low 
redshift, we assume a smooth intergalactic/intracluster 
medium (IGM/ICM) and gas-poor systems, i.e. no disk 
of c older (i.e. atom ic or molecular) material in the host 
(c.f. [Sutherland fc B ickncU 2007). We note that in sev- 
eral exploratory runs, additional atmospheric complexity 
is overlaid onto a smooth atmosphere with no change in 
the gravitational potential $, i.e. not initially in hydro- 
static equilibrium. 

The jets in our simulations are hypersonic (^lOOc^) 
light (pjct ~ O.Olpo) flows injected in pressure equilib- 
rium with the ambient material from back-to-back circu- 
lar footpoints on the small inner boundary sphere at the 
origin. The kinetic luminosity of the jets is given by 



(11) 



where ^jct is the area of the footpoint at rinnor, and 
ignoring higher order contributions from the gravita- 
tional energy or thermal flux. To tune Lkin we pri- 
marily vary ujet because (i) Lkin is most sensitive to 
changes in tijet, (ii) the jets must be "light" to ensure 
that the Kelvin-Helmholtz instability growth rate at the 
boundary of the cocoo n is approximately th e same as 
in the relativistic case (jRevnolds et al.l[2002t ). and (iii) 
jets are highly collimated. The maximum jet width is 
constraine d by high-resolution X-ra y observations of jet 
knots (e.g. iPerlman fc Wilsonll2005[ ) and very long base- 
line i nterferometry observa tions of transverse structures 
(e.g. IGabuzda et al.l [20041 ). However, we note that in 
Section 3.2 we vary pjet and Ajct in a limited range. 

We insist that the jet cover enough grid zones (~ 
30) to resolve transverse structures such as the oblique 
shocks which collimate the jet. The jet is rapidly pre- 
cessed around a small angle (a < 0jct) at 207r rad 
(code time) to break the axisymmetry of the sim- 
ulation setup and simulate helical instabilities, thus 
spreading the thrust o ut over a larger working sur- 
face fc.f. iVernaleo fc Re ynolds 2006: He inz et al. 120061: 
[Sutherland k. Bicknellll2007tiO'Neill fc Jonesll2010l) . The 
location of the jet on the inflow sphere is thus given by 
(6', 0) — (a, rijct • i)- Although the angular quantities are 
free parameters, they are constrained by the observed 
coUimation of jets. Typical values are 0jet = 7r/7.5 and 
a = 7r/60 for nnnor = 0.05. In agreement with other 
work, these jets develop a cylindrical core of fast-moving 
material sheathed in a slower concentric shell continu- 
ous at the boundaries with the surrounding material and 
core. 

We use a grid with zones spaced according to a geo- 
metric series in r and 9 (256 bins each) and uniformly in 
(64 bins). This grid resolves the internal jet structure 
and important processes near the injection footpoints. 
We choose r € [0.05, 5.0] and S [0.01, 3.13] to avoid the 
polar singularity and begin with reasonably sized grid 
zones. The 9 grid is broken into two symmetric 128- 
bin components with the smallest zones near the poles 
where the jets are injected. The adequacy of our grid 
is demonstrated by a run with double the resolution in 
each direction which produces similar internal jet struc- 
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ture. Mixing is not substantially different at the lobe 
boundaries. 

We use periodic boundary conditions in the (j) direction 
and reflecting boundary conditions in 6. The outer r 
boundary at r = 5.0 has outflow conditions (material 
leaving the grid); likewise, outflow conditions exist for 
all finncr zoues BxcBpt wherB the jets are injected. The 
outer boundary is far from the jet activity so negligible 
material leaves the grid there. We discuss the importance 
of the inner boundary sphere in Section 4. 

Code units are transformed to physical units by choos- 
ing ap propriate values for r p , Cg, and po- For in- 
stance, iRevnolds et"all pOOl defined tq = 100 kpc, 
Cs — 1000 km s~^ and = O.Olmn g cm^"^ for a 
rich cluster and rg = 10 kpc, Cg = 500 km s^^ and 
Pq — O.lmjj g cm^^ for a group or poor cluster. In 
the former scheme a code unit of time (derived from 
the crossing time) corresponds to 50 Myr and in the 
latter 10 Myr. In our runs, we vary rg and po but fix 
Cs\ a value of Cg ~ 500 km s~^ appears appropriate for 
XRGs based on temperatures derived from spect ral fit- 
ting (HodgeilKl^cEeFiniloiS Ean^^ A jet 
injected at lOOcg would then have a physical speed be- 
tween 0.17c — 0.35c (a Lorentz factor of 7 ^ 1.01 — 1.05). 
However, as noted bv IKomissarov fc Falld (|1996( ). com- 
paring classical and relativistic jet simulations requires 
careful matching of parameters, in particular the mass- 
energy density content of the jet. Hence, the jet veloci- 
ties chosen should not be taken directly as assumptions 
of true jet velocity. 

Our runs were parallelized and used variety of proces- 
sors. Many of the runs were conducted on quad-core 
Intel® Core^"" 2.4 and 2.83 GHz workstations. 

2.2. Strategy 

We now outline our guiding strategy to determine 
whether X-shaped sources can result from the interac- 
tion of radio galaxy lobes with their environment. 

First, we only model powerful FR II sources which pro- 
duce strong back- flows. Neither the physical origin of the 
FR I/II dichotomy nor the differences between jets in ra- 
dio loud and radio quiet sources are understood, but the 
higher luminosity FR lis exhibit the hot spots, well de- 
fined lobes, and bridges exploited by the backflow model. 

In accordance with the observations 

(ICapetti et all I2002t ISaripalh fc Subrahmanv an 2009; 
IHodges-Kluck et al.l 120101 ) . we inject jets along the 
major axis of the surrounding atmosphere. We first use 
unrealistic atmospheres with very favor able pressure 
gradi ents for wing production (as in ICapetti et al.l 
1200 3). then explore jet and atmosphere parameters to 
study the production and characteristics of XRGs. In 
particular, we explore the dependence of wing formation 
on jet width, density, and kinetic luminosity as a 
function of time, and initial atmosphere parameters 
(core radius, density, ellipticity, and position angle). 
Although some parameter combinations are degenerate, 
this is a very large parameter space because the Lkin(i) 
may include dead time and intermittency. Compared 
to 2D modeling, our 3D simulations eliminate the re- 
quirement of a xisymmetry which enhances the jet head 
advance speed ()Bodo et al.lllQQS l) and the coherence of 
the back-flows. The effect of turbulence is also more 



realistically explored in three dimensions. 

Motivated by the expectation that pressure gradients 
affect the wings fundamentally the same way in different 
systems, we start from the ansatz that jet and atmo- 
spheric parameters are orthogonal. In other words, the 
jet parameters may be tuned in some fiducial atmosphere 
and the atmosphere parameters may be tuned with some 
fiducial jet such that the behavior of an arbitrary jet in 
an arbitrary atmosphere can be inferred. On the basis 
of this method we will evaluate the factors key to wing 
prominence and the viability of the backflow models. 



2.3. Hydrodynamic Models 0/ Powerful Double-Lohed 
Radio Galaxies 

If wings are distortions to a generic double-lobed 
source, we expect our simulated radio galaxies to strongly 
resemble the double-lobed sources produced by earlier 
models. Hydrodynamic and MHD models of jets are 
commonly employed to study either the phenomenology 
of the jets and lobes (with relativistic and nonrelativistic 
fluids) or the energy deposited by the lobes into their en- 
vironment. Hence, a large body of work concerned with 
phenomenological models of these sources exists and has 
produced key insights into the life of a powerful radio 
source. 

These studies include work on early jet /lob e evo- 
lution in a relaxed a tmosphere (e.g. IKrausd I2Q05I : 
lAntonuccio-Delogu fc Silk ,2010,1 . the influence of the 
Pjet/po de nsity contrast and intern al Mach num- 
ber (e.g. ICarvalho fc O^eal [2002al: IKrausd 120031: 
iVernaleo &: Revnoldsl |2007[) , iet stability in nonrela- 
tivist ic and relativist i c con ditions fe. g. IRosen et al.l 
1991 iKe ppens et all I200I rO'Ncih fc Jones! 120101 : 



Mignonc et al. 201(T), the importance of intermittency 
([O'Neill fc Jones ,20101 ) or intrinsically spreading out 
the jet thrust to slow the advance of simulated jets (e.g. 
iHeinz et al.ll2006l) , the importance of the jet environment 
to morphology or energy deposition (e. g. Capctti et aL| 
200^; Carvalho & O'Dca''2002b'; 'Krause 2005; Zicr 2005; 
Vernaleo fc Reynolds .2007; .Simionescu et al.. , 2009 : 
Kawakatu et al.l 120091: iFalceta-Gon^alves et al.l I2OIQ : 



Gaibler et al 



2010D . and physics beyond MHD (e.g 



Saxton et al.l 20lW . Naturally, there is overlap between 



the phenomenological studies and those motivated by 
the challenge to produce jets which prevent a cooling 
catastrophe in the ICM of a host galaxy cluster by 
isotropizing energy distribution. The viability of radio 
galaxies as AGN feedback mechanisms is beyond the 
scope of this work, but jet lifeti me and the passive 
evolution of "dead" radio galaxies (R eynolds et al.|[2002D 
place important constraints on wing prominence. 

Based on this work, we understand a powerful double- 
lobed source to have three distinct phases of evolution: 
(i) ignition, in which a spheroidal cocoon of relativistic 
plasma is formed around the nascent jets, (ii) the active 
phase, wherein the jet produces the cigar-shaped lobes 
associated with FR II radio galaxies, and (iii) the pas- 
sive phase where lobes evolve after the jet is turned off. 
Below we outline these phases for a preliminary hydrody- 
namical simulation of a fast {vjot = lOOCg exp(— 3t)), light 
jet with large Lkin (a "FR H" source) to provide a frame- 
work for discussing the evolution of winged sources (left- 
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Figure 2. Comparison of a normal (left panels) and winged (right panels) radio source using the same jet {vjat ~ 100exp[— 3t]) 
in atmospheres differentiated only by e. Note the size of the overpressured cocoon relative to the size of the atmosphere in the 
top panels. At a code time of 1.0 (100 Myr) the jet was restarted (ujct = 100exp[— 3(t — 1.0)]). The colorbars are the same in 
each image; the integrated radio is brighter for the normal source because the winged source has the same amount of material 
spread out over more volume. 

acquire a circulatory motion and flow along the inner 
edge of the bow shock (Figure |3l top panel). By the 
time these flows reach the midplane between the two 
jets, they have velocity vectors pointing radially inwards 
and do not collide with back-flow s from the counter-jet 
(jAntonuccio-Delogu fc SillJ l2010( ). In our models, the 
back-flows are also prevented from interacting with their 
counterparts from the counter-jet because our jets are in- 
jected into pristine atmospheres such that at very early 
times, each jet inflates its own cocoon. The two cocoons 
are both bounded by a bow shock, and the two shocks 
meet at the midplane between the jets and form an inter- 
stice of dense material that prevents early mixing (Fig- 
ure [Sj top panel). Whether such an interstice is present 
in real sources is not clear (we do not resolve the very 
earliest jet stages such as the "fio od-and-channel" phase 
seen in a clumpy, warm disk in [Sutherland fc Bicknelll 
l2007f ). but the interstice is ablated and the cocoon is 



hand panels of Figure 12. 2|) . We note that iBelan et al.l 
(|2011t ) have recently found good agreement between the 
structures observed in hydrodynamic simulations of hy- 
personic jets and laboratory experiment. 

Ignition — The jet is injected in pressure equilibrium 
with the initially relaxed surrounding medium. The flow 
quickly forms and drives a spheroidal bow shock into 
the surrounding medium (Figure 12. 2[ top left panel). 
At the location where the jet impacts the shock, back- 
flows develop which fill the space evacuated by the bow 
shock with light, hot spent jet plasma. This plasma 
forms a cocoon which sheaths the jet. At very early 
times the jet expands laterally, since it is unconfined by 
the cocoon plasma (in part due to the initially conical 
shape of the jet) and the exp ansion of the nascent radio 
galaxy is nearly self-similar (jCarvalho fc O'Deal l2002al : 
[Krausai2005i ). However, as a result, the early back- flows 
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unified by the active phase. 

As the cocoon pressure builds due to backflow con- 
fined by the bow shocks, the jet becomes azimuthally 
confined and a recoUimation shock appears near the in- 
jection point; by the end of the ignition phase, the cocoon 
is strongly overpressured relative to its environment. At 
the same time, the terminal shock becomes increasingly 
distinct from the bow shock, and the jet head takes on 
the characteristic double-pronged appearance of a 3D hy- 
drodynamical jet. While the back-fiows still follow the 
bow shocks during the late ignition phase, they are in- 
creasingly straight. We note that this early evolution is 
seen in all our runs and does not depend much on atmo- 
sphere or jet parameters, nor on the boundary conditions 
(as long as the core radius of the atmosphere vq is much 
larger than nnner)- 

Active Phase — When the jet head overtakes and pierces 
the initial bow shock s (driving a b ow shock contiguous 
with the initial burst: lKrausell2005l ). the active phase be- 
gins and a classical radio galaxy develops, with a cigar- 
shaped cocoon, hot spots, bow shocks (middle-left panels 
of Figure [221) , ^-^^d straight back-fiows (Figure [3l bottom 
panel). At this point, the cocoon ceases to be substan- 
tially overpressured and the cocoon's lateral expansion 
falls close to the sound speed of the ambient medium, 
lagging behind the bow shock (the cocoon expands due 
to buoyancy, Kelvin- Helmholtz, and Ray leigh- Taylor in- 
stabilities). This weak shock (wshock '~ 1.5cs) sweeps up a 
large amount of material as it expands and becomes bet- 
ter described as a strong sound wave (Figure [121) ■ Mean- 
while, the jet develops internal oblique collimating shocks 
along its length, retaining a high velocity out to the ter- 
minal shock at the jet head (i.e. the angular size of the jet 
decreases with increasing radius) where the cocoon be- 
comes momen tum-driven. Less powerful j ets do not pro- 
duce cocoons ()Vernaleo fc Revnoldsll2007l) . but we focus 
on the cocoon-bounded case (the coco on pressure also de- 
pends on the internal Mach number; iCarvalho fc O'Deal 
[2002a). 

Passive Phase — Once the jets are turned off, the 
back-fiows cease and the "dead" radio galaxy rises 
buoyantly in the atmosphere, lifting large amounts of 
material to great heights and mixing w ith the sur- 
rounding medium ([Reynolds et al.l l2002| ). Although 
the radio galaxy is relatively efficient at depositing 
its energy irreversibly in the surrounding medium, 
it is difficult to isotropize this energy deposition on 
short enoug h timescales to avoi d a cooling catas- 
trophes (e .g. 'Reynolds et al.' '2002'; 'Omma et al.' '2004; 
Omma fc Binncv 2004; Vernalco & Reynolds 2006, 2003 
O'NciU et al. 2005; DeYound l2010t lO'Neih fc Jon^ 
"OlOc ,Ostriker et al.. ,20"Tol ). The cocoon separates into 
bubbles which pinch off along the direction of the 
jets and rise in opposite directions; the evolution of 
these bubbles may be quite complex (e.g. Begclman 
200 1"; 'Churazov et al." 2001'; 'Ruszkowski ct al.^ 2007; 
Dong fc Stone 2 Q09: O'NciU ct al. 2009; Braithwaiia 
2010HPope et al.ll20'l0h . This very late stage evolution 
is unlikely to be important for winged galaxy evolu- 
tion, as winged gala xies appear to be related to sources 
with strong bridges ([Leahy fc Willianislll984[ ). and is not 
shown in Figure 12.21 



3. RESULTS 

With the scheme described in Section 2, we produce 
winged radio galaxies by the deflection of backflow into 
channels perpendicular to the jets. In our models, wings 
are produced in strongly asymmetric atmospheres when 
the jet is driven near the major axis; the wings are 
produced along the minor axis as cocoon material es- 
capes the central regions of the atmosphere and evolve 
buoyantly, coUimated by the surrounding stratified at- 
mosphere which promotes expansion along the steepest 
pressure gradient. Hence, these wing s ar e similar to 
those produced bv lCapetti et al.l ([20021 ) and [Zanni et al.l 
([2003| ). Without introducing additional complexity to 
the atmosphere, hydrodynamic wings are formed in two 
stages. First, during the ignition phase (Section 2.3), 
the overpressured cocoon expands faster along the mi- 
nor axis due to the steeper pressure gradient and forms 
small channels (proto-wings) perpendicular to the jet 
into which back-flowing plasma flows. The proto-wings 
produced in this way can account for 20-40% of wing 
length at the end of the active phase depending on the 
gradient. Second, during the active phase, the wings 
rise buoyantly, fed by turbulent flows near the midplane. 
Although the back-fiows near the jet heads are initially 
laminar and supersonic (relative to the lobe material), 
they quickly decelerate and do not enter the wings as 
coherent fiows; the wings expand subsonically. 

During the active phase the wings evolve almost inde- 
pendently of the cocoon. Hence, once wings have devel- 
oped, their length depends only on the properties of the 
native atmosphere whereas the length of the jet-fed pri- 
mary lobes is dominated by the properties of the jet (in 
particular the kinetic luminosity as a function of time) for 
powerful jets. Therefore, decaying jets produce promi- 
nent wings. 

In this section, we explore these ideas in detail, first 
comparing the life of a winged source to a canonical 
double-lobed one (Figure 12. 2p and then looking at the 
dependence of wing prominence and longevity on vari- 
ous tunable parameters in our models. 

3.1. Evolution of a Winged Source 

Winged galaxies experience the same life stages out- 
lined in Section 2.3. We describe how wings fit into this 
process in detail below, referencing our standard (ex- 
tremal) atmosphere (Run STANDARD, Tables [1] and [SJ. 
While unrealistic, this atmosphere provides important 
insights into the backfi ow model and can be directly 
compared to prior work ([Capetti et al.[[2002[ : [Zanni et al.l 
[2003) . 

Ignition — During the ignition phase of a winged source, 
the anisotropic expansion of the overpressured cocoon in 
an aspherical atmosphere produces channels which will 
later become the wings (the proto-wings). This pressure- 
driven expansion is supersonic but brief, since the cocoon 
quickly reaches pressure equilibrium. Although these 
channels need not be produced by an overpressured co- 
coon, some channels must exist for wings to form. 

The degree of cocoon expansion in a given direction 
depends on the pressure gradient experienced in that di- 
rection, and hence on the atmospheric parameters. In 
particular, the eccentricity of the atmosphere effectively 
changes the core radius tq seen by the cocoon in different 
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Figure 3. We show here and Vy velocity vectors over- 
laid on density slices (taken ai (f) — 0) at different times for 
the simulation TX_E75_VE3_B5.0 (Table 3). We have cho- 
sen (arbitrarily) Cs = 500 km s~^ and ro = 50 kpc (corre- 
sponding to a code ro = 2.0) to represent a large elliptical 
galaxy. White vectors point radially outwards and cyan vec- 
tors point radially inwards (i.e. the back-flows). Vectors are 
only shown when the magnitude of v exceeds the sound speed 
of the lobe material (cgjobo = lOcs). Note that early on, the 
back-flows follow the contact discontinuity and are directed 
towards the inner boundar y, but later become straight as in 
lAntonuccio-Deloeu fc Sili3 "|2010n . 



directions. For example, along the z-axis (rx; the jet is 
driven along z) Equation 4 gives: 



dp 
dx 



3 cl 3/2 

"2 7^°Vff 



(12) 



where ro,eff = To{^ — eniax)/(l — f-x) as above. Since the 
pressure gradient is steeper for smaller ro,cff, higher el- 
lipticity along other axes promotes wing expansion: for 
tx = 0.0 and Cz = Cmax = 0.75, ro,cff = 0.25ro along the 
a;-axis. Conversely, if ex = Emax = 0.75 and Cz = 0.0 (the 
jet is pointed along the minor axis), ro^off = ^"0 and wings 
are suppressed. Of course, wing expansion depends on 
the actual pressure gradients rather than just the ratio 
along different axes; the base values rp and po determine 
whether wings can form (i.e. a highly elliptical atmo- 
sphere can have shallow pressure gradients along the mi- 
nor axis if it is very large) . Note that at very small radii 
f ^ ''o.eff (i-e. during the ignition phase), the pressure 
gradient steepens linearly with increasing r and at large 
radii (during the active phase, see below) the pressure 
gradient becomes shallower as r~^/^. Hence, the ignition 
phase is the time at which atmospheric asymmetry has 
the strongest effect on the ultimate morphology of the 
source. 

Active Phase — Once the cocoon has reached pressure 
equilibrium, the proto-wings are no longer pressure- 
driven. As the back- flows become straight, they merge 
and fill the channels (Figure |3l bottom panel) to form 
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Figure 4. Positions and velocities of the leading edges of 
the bow shock, wings, and lobes shown for (a) the standard 
simulation and (b) a triaxial simulation. The dashed line in 
the top panel shows the position of a point moving radially 
outwards at the sound speed. Note the different a;-axes; at t = 
1.0 in the left panels, a new jet is reinjected. Negative growth 
for the lobes indicates collapse during the passive phase. 



structures more closely resembling observed wings. How- 
ever, the flows into the wings are turbulent and transonic 
or subsonic (relative to the internal lobe sound speed 
Cs.iobc = lOcs). In other words, wings are not driven 
during the active phase. Rather, because they are struc- 
tures filled with light fluid, they rise buoyantly in the 
atmosphere, collimated by the stratified elliptical atmo- 
sphere. 

Since the wings rise buoyantly, their growth rate is 
subsonic. To see this, it is instructive to look at the 
simplified case of a spherical bubble (of fixed radius) ris- 
ing buoyantly in a dense fluid. The (terminal) buoyant 
velocity at a given height for such a bubble is 



^buoy 



2.9 Vb 

CdA^ 



16 rb 



37 r^. 



cff - 



(r/rccff)^ 



(13) 



where Vb, ^b, and rb are the volume, cross-sectional area, 
and radius of the bubble respectively, Cd = 0.75 is the 
"drag" coefficient, and we have used the background dark 
matter potential <I> to compute g in terms of r and effec- 
tive core radius ro^eff = ''o(l ~ emax)/(l — e)- Naturally, 
this follows the pressure gradient (Equation 12), so when 
the bubble is at very small radii its velocity increases as 
r^^'^, whereas outside ro^off it decreases as r~^/^. For 
typical values in our simulations, the peak value of Ubuoy 
is approximately c^. Hence, in a highly elliptical atmo- 
sphere, a bubble rising along the minor axis will rise sub- 
sonically with decreasing velocity for most of its lifetime. 

Our wings do not fit this simplified case because they 
are large relative to ro^off, variable in size, aspherical, 
and connected to the cocoon. Nonetheless, their expan- 
sion is subsonic for the same reason: the wings exceed 
fQ.cff early on when their growth is pressure-driven rather 
than buoyant, and i^buoy monotonically declines there- 
after. This transition from supersonic to subsonic wing 
expansion is plainly seen in Figure HI 

Therefore, wings in our models may not intrinsically 
exceed the length of the primary lobes unless the primary 
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lobes advance at some average speed Vh < Cg . Since pow- 
erful jets are required to produce the proto-wings in an 
overpressured cocoon phase, prominent wings require de- 
caying jets, intermittent jets, or jets which deposit their 
thrust over an increasingly large area with time. The 
profile of the fiducial run is shown along with its wing 
prominence in Figure [S] and compared to a triaxial at- 
mosphere in Figure m Note that in the fiducial run (left 
panel of Figure]!]) the wings eventually overtake the lobes 
(in part due to the collapse of the lobes when the jet is 
very weak) , but clearly move into a subsonic regime early 
and remain there. The bow shock remains mildly super- 
sonic through most of both simulations. In Figure [S] we 
show wings produced by decaying jets during the active 
phase in several environments (the fiducial run is shown 
in the top row). 

Passive Phase — Once the jets are turned off, the cocoon 
disintegrates as the lobes either rise buoyantly as bubbles 
or collapse under the relaxing atmosphere. If the radio 
galaxy is within a few tq in a dense atmosphere, the fall- 
back of displaced material shreds the lobes into small 
bubbles. In either case, the wings pinch off and rise. 
These bubbles do not survive long, but their behavior 
might be substa ntially altered in the presence of mag- 
netic fields (e.g. lO'Neill et afl 120051 : iBraithwaitd IMl 
iPope et aLll2010D . Because observed winged sources have 
strong bridges, we expect them to be in the active phase. 

Remjection — A powerful jet reignited during the pas- 
sive phase before the cocoon has disintegrated may sig- 
nificantly enhance the wings. If the old jet channel is 
somewhat broken up, the reinjected jet forms a new ter- 
minal shock and bow shocks inside the old cocoon. These 
new shocks do not form spheroidal structures (as in a re- 
laxed atmosphere) but instead produce strong, straight 
back-fiows near the midplane. Hence, the wings receive 
a large infiux of fresh supersonic plasma directly after 
reinjection (Figure 12. 2[ bottom panels) . This brightens 
the wings substantially and reinforces their structures. 

There is a relatively narrow window of time where this 
process is effective. If the reinjection occurs while the jet 
channels are largely intact (i.e. during the active phase), 
the jet simply follows these channels. On the other hand, 
if the reinjection occurs when the wings have already sep- 
arated from the lobes as individual bubbles, the new jet 
cannot feed them. Even if the reinjection occurs at the 
"right" time, the efficacy of the restarted jet at promot- 
ing wings is short-lived. However, this mechanism can 
produce wings that are intrinsically longer than the jet- 
driven lobes (Figure 3, bottom panels). If the reinjected 
jet decays, this extreme axial ratio can be maintained for 
most of the lifetime of the restarted radio galaxy and a 
bona-fide X-shaped source results, although this source 
would only have the "FR H" primary lobes for a short 
period of time (~ 5 Myr for a source 40 kpc across). 

3.2. Ecology of Winged Sources 

In Section 3.1, we describe the life stages of a winged 
radio source. We now describe the dependence of wing 
prominence on our tunable parameters. Following the 
strategy laid out in Section 2.2, we begin with the fidu- 
cial simulation STANDARD which is characterized by an 
atmosphere with a large ellipticity e and a jet with a 
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Figure 5. Plot of W / Lmsx (crosses) in timestep intervals of 
5t = 0.05 for the standard simulation with a reinjection at 
t — 1.0. We have overplotted the velocity of the jet as the 
black line (right j/-axis). The Roman numerals indicate (I) the 
overpressured cocoon phase, (II) the active phase (powerful 
jet), (III) the active phase (weak jet), (IV) the passive phase 
(with cocoon collapse), (V) a second active phase due to the 
reinjected jet, and (VI) the second passive phase with further 
cocoon collapse. During cocoon collapse, the source does not 
resemble a winged source (instead it is a "dead" radio galaxy). 
A source appears winged during the active phase when the 
jet is most like a weak FR II or strong FR I. 




Figure 6. Examples of different intrinsic wing morphology. 
Left: Initial relaxed atmosphere density. Center: Winged 
sources during the late active phase. Right: False synchrotron 

maps {ji, oc pI p'^^ assuming equipartition) of these sources. 



Hydrodynamic Models of Radio Galaxy Morphology 



9 



velocity profile ujet = lOOcg exp[— 3t] (Figures 12.21 and 
|4]). Taking the jet from the fiducial run, we individually 
vary parameters in the atmosphere in order to see their 
influence on wings (Figure [71 Table [IJ , and then vary 
jet parameters in the fiducial atmosphere (Figure |H1 Ta- 
ble [5]) . We then synthesize the information gleaned from 
these single-parameter curves to attempt to find atmo- 
spheres and jets which can explain real winged sources, 
including more complex behavior as well (Table [3]) . All 
of the runs proceed to a code time of 1.0 or more, where 
most have entered the passive phase. Reinjections, where 
used, occur near the beginning of the passive phase. 

To quantify wing prominence we use the wing-to-lobe 
axial ratio {W/L). This quantity is an imperfect measure 
because W/ L intrinsically varies with azimuthal angle </> 
and time, and observable sources are seen in projection 
so the observed W/L will differ from the intrinsic value. 
W/L is also essentially meaningless during ignition or 
the passive phase, but a source in one of these stages 
would not be classified as "winged." For the parameter 
exploration in Figures [7] and [8] (values in Tables [T] and 
[5]), we adopt as a fiducial value VF/Lo.s^ the maximum 
intrinsic W/L for any pair of azimuthal angles + tt] 
at a code time t — 0.5. For the standard jet (wjct = 
lOOCs exp[— 3t]), t — 0.5 represents the transition from 
a powerful to weak jet in the active phase, and by this 
time the wings have begun to grow buoyantly. Hence, 
a code time of t = 0.5 is a reasonable place to measure 
the infiuence of the atmosphere on wings. We also use 
in varying the jet parameters, noting that while 
t = 0.5 is no longer special, all of the runs in Table [5] are 
in the active phase at this time. W/Lq,^ is not predictive 
of wing length later in the same simulation, but is a good 
measure of relative wing prominence between simulations 
due to the subsonic growth of wings. 

From Figure [7] one can get a broad sense of the de- 
pendence of wing prominence on the size and shape of 
the atmosphere. It is immediately obvious that there 
is a strong dependence on the ellipticity e of the atmo- 
sphere, and that smaller, denser atmospheres are the 
most conducive to wing formation (although "small" 
could be physically quite large depending on the jet). 
It is also notable that relatively high values of W/L 
can be achieved hy t = 0.5; because we use a de- 
caying jet, the wings will only become more promi- 
nent during the later active phase. While these re- 
sults are not surprising (Equation 12), the particular 
form of the curves depends on both the jet and the 
atm osphere. As expected fro i n iCapetti et al.l (|2002[ ) 
and iSaripalli &: SubrahmanvaiJ ()2009l ). long wings re- 
quire jets co-aligned close to the major axis. 

The trends are less clear when varying jet parame- 
ters (Figure In the top panels, the jet velocity is 
allowed to vary freely without conserving integrated ki- 
netic luminosity between runs. If the jet is injected at 
constant velocity throughout the simulation (Figure [Sti), 
a faster jet is slightly better than a slower one at mak- 
ing wings. However, a decaying jet (Figure ^p) is bet- 
ter still, with pure exponential decay more effective at 
increasing W/Lq,^ than "Gaussian" jets of the form 
Vjet = voexpl—at"^]. In the bottom panels of Figure [51 
we hold the kinetic luminosity constant between runs 
by modifying vq. From Figure [5};, it is evident that in- 



Table 1 

Varying the Atmosphere using the Standard Jet 



Name Wings? ro po APA W/L 
(deg.) (f = 0.5) 

Standard Atmosphere 

STANDARD Y 1.0 3.0 0.0 0.75 0.0 0.62 



EUipticity 



SJ.E20 


N 


1.0 


3.0 


0.0 


0.20 


0.0 


0.42 


SJ.E30 


N 


1.0 


3.0 


0.0 


0.30 


0.0 


0.43 


SJ.E50 


N 


1.0 


3.0 


0.0 


0.50 


0.0 


0.49 


SJ.E60 


Y 


1.0 


3.0 


0.0 


0.60 


0.0 


0.54 



Core Radius 



SJ.R0.5 Y 0.5 3.0 0.0 0.75 0.0 0.73 

SJ.R0.75 Y 0.75 3.0 0.0 0.75 0.0 0.61 

SJ_R2.0 N 2.0 3.0 0.0 0.75 0.0 0.45 



Core Density 

SJ.D1.5 Y 1.0 1.5 0.0 0.75 0.0 0.52 

SJ.D2.0 Y 1.0 2.0 0.0 0.75 0.0 0.57 

SJ.D4.0 Y 1.0 4.0 0.0 0.75 0.0 0.65 

SJ.D5.0 Y 1.0 5.0 0.0 0.75 0.0 0.68 



APA 



SJ.PA5 Y 1.0 3.0 0.0 0.75 5.0 0.62 

SJ.PA10 Y 1.0 3.0 0.0 0.75 10.0 0.60 

SJ.PA15 Y 1.0 3.0 0.0 0.75 15.0 0.58 

SJ.PA20 Y 1.0 3.0 0.0 0.75 20.0 0.51 



Triaxial Atmospheres 

TXJ]75.1 Y 1.0 3.0 0.375 0.75 0.0 0.65 

TXJ]75.2 Y 1.0 3.0 0.50 0.75 0.0 0.72 

Note. — Runs in which the standard atmosphere was 
varied one parameter at a time, holding the standard jet 
(«jet = lOOcs exp[— 3t], a = 7r/35 rad and /3 — n/15 rad for 
dinner = 0.1) coustant. The STANDARD run uses the standard 
jet and atmosphere and hence is a data point in each cate- 
gory. The "wings" column denotes whether a run produced 
noticeable wings at any point during its active lifetime. W/L 
at t = 0.5, on the other hand, is a way to directly compare 
different runs. At t = 0.5, wjet ~ 20cs, i.e. twice the lobe 
material sound speed and a transition point during the ac- 
tive phase between a powerful and weak jet. The time of this 
transition depends on the velocity profile of the jet; t = 0.5 is 
only correct for the standard jet. APA is the angular distance 
between the jet and the major axis. 

creasing the width of the jet (at the cost of a slower 
jet) is important. Changing the density of the jet mate- 
rial has an extreme effect on M^/Lq.s (Figure [8ll) ; above 
Pjot ~ 0.01, the jet ceases to be "light" relative to the 
backgrou nd and no longer form s lobes resembling radio 
galaxies (jRevnolds et al.ll2002D . At very low densities, 
mixing becomes very efficient at disrupting the radio 
galaxy. 

We now discuss these results in detail, beginning with 
the properties of the fiducial run. 

Fiducial Run— Our fiducial run (STANDARD, Tables[I]and 
[5]) is the combination of our fiducial atmosphere (rg = 
1.0, Po — 3.0, €x ^ £y ^ 0.0, = 0.75) with our fiducial 
jet (wjet = lOOcs exp[— 3i], 6'jct = 7r/15, a = tt/35, aligned 
along the z-eods). Like most of our runs, this simulation 
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Figure 7. These plots show the dependence of axial ratio 
W/L on (a) the ellipticity of the atmosphere e, (b) the core 
density po, (c) the core radius ro, and (d) the angle between 
the jet and the major axis APA. Each plot represents varying 
one parameter while holding the other standard atmosphere 
parameters (e = 0.75, po ~ 3.0, ro — 1.0, and APA = 0) 
steady and using the standard jet. In panel (a) the square and 
triangle represent runs with different ex for = 0.75, demon- 
strating the effect of triaxiality on W/L (Table [!}. W/L is 
measured a,t t — 0.5 for all cases; for the standard jet, this 
represents the transition during the active phase from a pow- 
erful to a weak jet (this time and the "wing threshold" vary 
with choice of jet). Values taken from Table [1] 



proceeds to a code time of 1.0, at which point the radio 
galaxy is in the early passive phase. Unsurprisingly, this 
run produces some of the most prominent wings of our 
suite (Figure [6]); the shallow pressure gradient along the 
2-axis and declining velocity profile combine to produce 
stalled lobes which begin to collapse at the end of the 
active phase. At the same time, wings quickly escape the 
core region, and hy t = 0.4 are easily recognizable (right- 
hand panels of Figure [2?2t . After reaching W/L = 0.83 
at t = 0.75, the cocoon collapse leaves the wings as the 
most notable features (Figure [5]) . 

At i = 1.0, we re-inject a jet with Vjct = lOOcs exp[3(i — 
1.0)] (Figure [51 and bottom panel of Figure Since 
the jet is expanding into tenuous material from the old 
cocoon, the expansion quickly re-establishes a cocoon. 
However, the reinjection also drives a bow shock inside 
the collapsing cocoon, allowing initially strong back-flows 
to feed the wings directly. Hence, between t = 1.10 and 
t = 1.35, an X-shaped radio galaxy is apparent {W/L ^ 
1.2). After t = 1.35, the jet again weakens to the point 
where collapse begins. If the jet is instead reinjected at 
the inception of the passive phase at t 0.75, the end 
result is very similar {W/L ~ 1.2 at t = 0.95). 

The atmosphere i n this run is si n iilar t o the 2D simu- 
lations pr esented in iCapetti et al.l (|2002D and 3D simu- 
lations in iZanni et al.l (|2003[ ). Striking asymmetries de- 
velop despite the axisymmetric atmosphere because of 
turbulent mixing between different slices in (j) and be- 
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Figure 8. These plots show the dependence of axial ratio 
W/L on (a) jet velocity for a constant-velocity jet, (b) e- 
folding time for a decaying jet (diamonds represent jets of the 
form Vjct = 100exp[— at] and triangles Vjct = 100 exp[— at^]), 
(c) area of the jet nozzle, and (d) density of the jet material. 
All runs are conducted in the standard atmosphere (Table [1]). 
In (a) and (b), the jet kinetic luminosity is not conserved 
between runs, whereas in (c) and (d) we adjust the initial 
velocity no to conserve Lkin- As in Figure [T] we show W/L as 
measured at a code time t = 0.5; all sources are in the active 
phase but not at the same place. 



cause the back-flows form three-dimensional structures 
within the lobes. If these flows are slightly misaligned on 
opposite sides of the midplane, asymmetries develop. 

Atmosphere — In our exploration of different atmosphere 
properties (Figure [T] Table [l]), we vary e, po: ''Oi and 
APA individually while keeping the other atmosphere 
parameters steady and using the standard jet {vjct = 
lOOcs exp[— 3t], pjct = 0.01, a = n/35, and (3 = 7r/15). 
Except when we vary APA, the jet is coaligned with the 
major axis. 

It is clear from Figure [7] that wing prominence depends 
strongly on the core radius tq, core density po, and el- 
lipticity e. Because the lobes and wings evolve almost 
independently, these dependencies can be understood by 
breaking the axial ratio W/L into its constituent parts: 
the wing length W and lobe length L. For the purposes 
of this discussion, W and L are both taken to be mea- 
sured at t = 0.5 (Table [J). 

We find that L varies only slightly as a function of tq 
(less than 10%) but does peak between rg — 0.5 and 
ro = 2.0. At very low tq, the lobes are not confined by 
the atmosphere and spread out laterally, causing more 
turbulent mixing and dissipating the thrust of the jet. 
At very high tq, the lobes are completely confined by 
the atmosphere and the density gradient is so shallow 
that as the jet weakens it is increasingly resisted by the 
atmosphere. The importance of this effect on the lobes 
is more pronounced in a comparison between L and po 
where the lobe length dechnes by 50% between po = 1.5 
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and po = 5.0. 

On the other hand, the wing length W depends much 
more strongly on the core radius tq (declining precipi- 
tously with increasing rp) and only weakly on the core 
density. This behavior occurs because the wings re- 
spond primarily to the pressure gradient along the mi- 
nor axis (Equation 12) which depends on and po as 
VP oc porQ^. Thus, a large W/L requires small, dense 
atmospheres. 

In this context, the importance of ellipticity is clear: 
e determines the ratio of core radii along the major and 
minor axes. For high e, the lobes expand into material 
which does not vary much in density during the active 
phase whereas the wings quic kly escape the c ore re gion. 
This behavior is unsurprising; iCapetti et al.l ()2002 f) and 
[Kraft ct al. (2005) make the same basic argument for 
wings produced when the jet is coaligned with the major 
axis of the atmosphere. However, the interplay between 
the parameters rules out predicting the wing length from 
first principles. For instance, one might imagine that we 
could produce wings in an atmosphere with e ~ 0.3 (Fig- 
ure [7]) by making vq tiny and po large. Although this 
does increase the wing length, the radio galaxy expands 
beyond the dense part of the atmosphere during the ig- 
nition phase so the jet does not face as much resistance 
as might be expected. 

It is possible to produce wings in atmospheres with 
smaller e by making the atmospheres triaxia l (using pro- 
late e llipsoids) instead of axisymmetric. IZanni et al.l 
(|2003D expect triaxial ati nospheres (as oppose d to the ax- 
isymmetric simulations of lCapetti et al.lf2002D to produce 
longer wings because the backfiow is collimated along a 
single minor axis as opposed to forming a torus at the 
midplane. We find that while a triaxial atmosphere does 
collimate the wings, it does not by itself make them sub- 
stantially longer relative to an axisymmetric atmosphere 
with the same maximum e because the backfiow in our 
simulations is a compressible fluid and wing growth is 
driven by buoyancy (the pressure gradient is the same 
along the minor axis). Rather, triaxial atmospheres pro- 
duce better-defined proto-wing channels during the igni- 
tion phase. These channels are then reinforced during 
the active phase. In other words, triaxiality increases 
W/Lo,5 at smaller values of e because the ignition co- 
coon is less axisymmetric. 

Finally, we test the sensitivity of wing formation to the 
degree of alignment between the jet and major axis APA 
(Figure [TH)- Long wings are difficult to produce when 
the jet is misaligned with the major axis, and the wings 
produced differ in character. This is because the ignition 
cocoon produces proto- wings perpendicular to the jet in- 
stead of along the atmosphere's minor axis. Therefore, 
the longer wings which develop during the active phase 
do not benefit as much from the initial supersonic ex- 
pansion. Wings also become increasingly associated with 
the lobe on their side of the major axis as APA increases 
because the lobes are bent by the atmosphere and the 
back-fiows no longer make it all the way to the midplane 
before flowing into the wings. Hence, the lobe- wing pairs 
are mirror-symmetric about the midplane. The angle be- 
tween the lobe and the wing in each pair is also largely 
determined by APA. 



Table 2 

Varying the Jet in the Standard Atmosphere 



Name Wings? pj^t Vj^tit)/cs a /3 W/L 

[t = 0.5) 

Standard Jet 



STANDARD Y 0.01 100exp[-3t] 7r/35 n/15 0.62 



Velocity Profiles 



SA.V30 


N 


0.01 


30 




■it/35 


7r/15 


0.39 


SA.V50 


Y 


0.01 


50 




it/35 


7r/15 


0.43 


SA.V80 


Y 


0.01 


80 




7r/35 


7r/15 


0.44 


SA.VlOO 


Y 


0.01 


100 




7r/35 


7r/15 


0.47 


SA.VEl 


Y 


0.01 


100exp[- 


■t] 


7r/35 


7r/15 


0.48 


SA.VE6 


Y 


0.01 


100exp[- 


■6t] 


7r/35 


7r/15 


0.67 


SA.VG3 


Y 


0.01 


lOOexpi- 


■3*2] 


7r/35 


7r/15 


0.47 


SA.VG6 


Y 


0.01 


100exp[- 


■6*2] 


7r/35 


7r/15 


0.51 


Jet Density 


SAJ30.005 


Y 


0.005 


100/exp[ 


-3t] 


7r/35 


7r/15 


0.73 


SAJ)0.05 


N 


0.05 


100/exp[ 


-3t] 


7r/35 


7r/15 


0.19 


SA_D0.10 


N 


0.10 


100/exp[ 


-3t] 


7r/35 


7r/15 


0.09 


Jet Width 


SA.B30 


Y 


0.01 


100/cxp[ 


-3i] 


7r/35 


7r/30 


0.51 


SA.B7.5 


Y 


0.01 


100/cxp[ 


~3t] 


7r/35 


7r/7.5 


0.76 


SA.B5.0 


Y 


0.01 


100/exp[ 


-3t] 


7r/35 


7r/5.0 


0.87 



Note. — The standard atmosphere is a relaxed, isother- 
mal /3-model with core radius ro — 1.0, core density po = 3.0, 
ellipticity e = 0.75, and with the jet oriented along the major 
axis. Few runs are needed to deduce the dependence of the 
morphology on jet parameters other than the kinetic lumi- 
nosity as a function of time, but in these cases we vary vo by 
some factor / such that Lkin(t) is the same as in the standard 
atmosphere; see text for caveats of this approach. 

Jet — In contrast to the atmosphere parameter explo- 
ration, determining the dependence of W/L0.5 on the 
character of the jet is difficult because of the potential 
for time dependence in the jet power. This time depen- 
dence, along with the small width of the jets, also makes 
it more difficult to test our findings observationally. We 
insist that the jets be light, hypersonic flows (Section 2) 
which produce sources which resemble double-lobed ra- 
dio galaxies, and within these constraints test the de- 
pendence of W/Lo,Q on jet power as a function of time, 
jet width, and density of jet material (Figure [8]). For 
each of the simulations in Figure |8] we use the standard 
atmosphere with ?'o = 1.0, po = 3.0, e = 0.75, and 
APA = 0.0°. Note that Figures [Hli and Figures [Hb do 
not conserve kinetic luminosity between runs, whereas 
Figures [8fc and Figures [8ji do. The exploration of jet 
velocity is conducted with jets with the standard width 
of /3 = 7r/15 (~ 0.01 code units wide) and density of 
Pjct = 0.01. 

Figures [5^ and [BJ^ demonstrate the importance of de- 
caying jets to long wings. In Figure [8^, we use jets 
with differing velocities but no time dependence, find- 
ing that none produce a large W/Lq ^. W/Lo,^ actually 
increases with increasing Vjct even though the faster jets 
also punch through the atmosphere more quickly. This is 
because the overpressured cocoon produced by a weaker 
jet during the ignition phase is less overpressured and 
thus produces smaller proto- wings. Below Vjet ~ 20cs, 
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the sources are not cocoon-bounded at all. Since the jets 
do not decay, W/L obviously decreases with time. 

Powerful jets which decay (Figure [Hd) are effective at 
producing wings because the proto-wings form during 
the ignition phase and the jet generates strong back-flows 
early on when the jet head is close to the midplane. Since 
the lobes grow increasingly slowly, the subsonically ex- 
panding wings keep pace with the active lobes more eas- 
ily. As is clear in Figure [8)d, the slower the decay, the 
smaller W/Lq,q (diamonds represent exponentially de- 
caying jets and triangles jets with wjet = exp[— at^]). 
Jets with increasing velocity do not produce long wings 
because a powerful jet is necessary at the inception of 
activity to produce proto-wings. 

In Figures IHt and Figures EJi we investigate the de- 
pendence of on jet width and density. We vary 
/3 and pjot while conserving Lkin by varying the ini- 
tial velocity of the jet wq, using the standard form of 
Vjct = Voexp[—3t]. We limit these runs to those with vq 
capable of producing an overpressured cocoon that seeds 
proto-wings. 

The nozzle width of the jet clearly has a strong effect 
on W/Lo,5 (Figure Ht). Increasing the jet width makes 
the ignition cocoon larger and more overpressured during 
the ignition phase, promoting expansion along the minor 
axis. The lobe width during the active phase increases 
with increasing /3, dissipating the jet thrust over a larger 
solid angle. Short, fat lobes attached to fat wings result. 
On the other hand, if /3 is tiny, the jet drills through 
the surrounding atmosphere quickly, and thin, long lobes 
result. It is not clear whether jet widths actually vary 
substantially between sources and what determines the 
width of the jet; all jets are very narrow. Thus, our 
results are more generally a statement that jets which 
dissipate their thrust over a wider area produce longer 
wings. 

Figure EJi shows W/Lo,^ as a function of pjot- The 
sharp decline with increasing pj^t is due to the thrust 
carried by the jet. Even at moderate velocities, denser 
material drives the lobes forward much faster than light 
material while at the same time driving weaker back- 
flows. Hence, wing formation is not favored. We note 
that when pjct > 0.05, the sound speed of the lobe mate- 
rial is not much greater than that of the ambient medium 
and the jet does not develop the usual KH instabilities. 
These jets are therefore not "light" as required to repro- 
duce realistic radio sources with non-relativistic hydro- 
dynamics (Section 2). However, these runs are shown 
along with the light ones to illustrate the importance of 
dissipating jet thrust to long wings. 

Synthesis Runs — We have used the insights gained from 
examining the dependence of wing prominence on various 
atmosphere and jet parameters to new runs with complex 
atmospheres and jets in an attempt to produce long, re- 
alistic wings. These simulations (listed in Table [3]) are 
not a systematic exploration of any phenomenon. We 
offer a few brief observations here. 

Triaxial atmospheres offer the best hope of making long 
wings at lower e (e.g. Figure |6l), but still cannot pro- 
duce substantial wings in our simulations for e < 0.45. 
Several runs from a large suite of such simulations are 
listed in Table |3] with the prefix TX. Generally, they fol- 
low the same trends as described in above; the impor- 



tance of both favorable atmosphere and jet parameters 
to wings is obvious in the simulation with ez = 0.45 
(TX_E45_D5.0_WIDE), where a wide, decaying jet in a 
small, dense atmosphere is required to produce wings 
comparable to that in the standard atmosphere at lower 
ellipticity. Embedded disks of hot material (i.e. a thick 
disk; runs DISK_VE1 and DISK_VE3), however, alleviate 
the problem, allowing prominent wings to grow in atmo- 
spheres with globally small e. These disks are somewhat 
denser than the larger ISM, so the ignition cocoon ef- 
fectively encounters a small, dense, highly elliptical at- 
mosphere during the ignition stage. The disks are blown 
apart by the blast wave from the ignition stage and might 
not be observable once the radio galaxy has turned on. 

In most of our models, the ignition stage occurs in 
a smooth, relaxed medium. In real galaxies, the ig- 
nition stage would occur in the galactic center where 
mergers, dynamical effects, and other phenomena associ- 
ated with AGN can significantly disturb the ISM. We do 
not attempt to model real galaxies, but find that adding 
Kolmogorov spectrum turbulence to our atmospheres or 
small bar-like perturbations to the underlying gravita- 
tional potential near the nucleus does not have a large 
effect on the ignition stage. As long as there is a suf- 
ficient amount of ISM distributed around the nucleus, 
an overpressured cocoon can form. On the other hand, 
in runs where we begin with a weak jet of Vjct = 20cs 
(twice the sound speed of the lobe material), and no co- 
coon is formed, wings do not form (the SLOW_START runs 
in Table [3]). The backflow model has difficulty producing 
wings when a weak AGN suddenly becomes more pow- 
erful because an overpressured cocoon only forms when 
the jet is completely confined by the atmosphere. These 
runs can, however, produce something akin to Z-shaped 
morphology by virtue of the lobes escaping the densest 
regions of the atmosphere before the powerful jet turns 
on. 

We have briefly investigated the potential for gas-rich 
"stellar shells" from minor mergers located periodically 
along the major axis to produce wings (motivated by 
IGopal-Krishna fc Wiital '20101. We use a simple model 
in which only the positive peaks of two-dimensional 
sine/sine waves are added to the density of an under- 
lying elliptical atmosphere; the center is cut out. The 
amplitudes are set at a maximum of po- The shells do 
not significantly impact the formation of wings emanat- 
ing from the center of the galaxy, although they do resist 
the jet, allowing for a hi gher W/ L (the SHELL run s in Ta - 
ble[31). In the model of iGopal-Krishna k Wii"tal (|2010D . 
wings are instead produced near the site of the shell, and 
we cannot reproduce this. Obviously, these runs are not 
a thorough exploration of the effect of stellar shells, es- 
pecially if they actually bend jets (in which case they are 
beyond the scope of our models) . 

Motivated by the importance of a reinjected jet (Sec- 
tion 3.1), we examine the role of rapid intermittency 
(runs INTERMITTENT in Table [S]). By this we mean jets 
which experience multiple outbursts during the run and 
which are essentially in an "on" state or an "off" state. 
Generally, we find that rapid intermittency has the effect 
of suppressing wing formation and results in much more 
regular cocoons (in terms of their projected morphol- 
ogy) than single outbursts. On the other hand, ragged 
cocoons are produced by long periods of dead time. In 
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Table 3 

Synthesis Runs 



Name 


ro 


PO 


ex 




Pjot 




Cs 


a 


/3 


W/ Lmax 


Interpolation 


WIDE.VEl 


1 


.5 


4.0 


0.0 


0.75 


0.01 


100 


exp 


-t] 


7r/90 


7r/7.5 


0.58 


WIDE.VE3 


1 


.5 


4.0 


0.0 


0.75 


0.01 


100 


exp 


-3t] 


7r/90 


7r/7.5 


0.73 


TX.E75.VE1 


1 


.0 


3.0 


0.375 


0.75 


0.01 


100 


exp 


~t] 


7r/90 


7r/7.5 


0.63 


TX.E75.VE1.BIG 


1 


.5 


4.0 


0.375 


0.75 


0.01 


100 


exp 


~t] 


7r/90 


7r/7.5 


0.68 


(restart at t = 1.0) 


















0.78 


TX.E75.VE1.B7.5 


0. 


75 


3.0 


0.375 


0.75 


0.01 


100 


exp 


't] 


7r/90 


7r/7.5 


1.02 


TX.E75.VE3.B5.0 


1 


.5 


4.0 


0.375 


0.75 


0.01 


100 


exp 


-3t] 


7r/90 


7r/5 


0.82 


TX_E60 


0. 


75 


3.0 


0.30 


0.60 


0.01 


100 


exp 


-t] 


7r/90 


7r/7.5 


0.44 


TX.E60.WIDE 


0. 


75 


3.0 


0.30 


0.60 


0.01 


100 


exp 


-t] 


7r/90 


7r/5 


0.67 


(restart at t = 1.0) 


















0.97 


TX_E60.D5.0 


0. 


75 


5.0 


0.30 


0.60 


0.01 


100 


exp 


-t] 


7r/90 


tt/7.5 


0.71 


TX_E45_D5.0_WIDE 


0. 


75 


5.0 


0.225 


0.45 


0.01 


100 


exp 


-3t] 


7r/90 


tt/5 


0.58 



Other 



DISK.VEl 


2.0 


3.0 


0.25 


0.50 


0.01 


100 exp 


-t] 


7r/90 


7r/7.5 


0.60 




0.25 


8.0 


0.0 


0.90 








DISK.VE3 


2.0 


3.0 


0.25 


0.50 


0.01 


100 exp 


-3t] 


7r/90 


tt/7.5 


0.83 




0.25 


8.0 


0.0 


0.90 








TURBULENT'' 


1.0 


3.0 


0.0 


0.75 


0.01 


100 exp 


-3t] 


7r/90 


7r/7.5 


0.64 


SHELL.VEl'' 


1.0 


3.0 


0.0 


0.75 


0.01 


100 exp 


-t] 


7r/90 


7r/7.5 


0.53 


SHELL.VE3'= 


1.0 


3.0 


0.0 


0.75 


0.01 


100 exp 


-3t] 


7r/90 


7r/7.5 


0.71 


INTERMITTENT 


2.0 


2.0 


0.0 


0.75 


0.01 


100 sin2 


[50t] 


7r/90 


7r/7.5 


0.44 


INTERMITTENT2 


1.0 


3.0 


0.0 


0.75 


0.01 


100 sin^ 


[6i] 


7r/90 


7r/7.5 


0.47 


SLOW.START.l'^ 


1.0 


1.5 


0.0 


0.75 


0.01 


20 to 


7r/90 


7r/7.5 


0.42 














100 exp 


-t] 




SLOW.START.2'^ 


1.0 


1.5 


0.0 


0.75 


0.01 


20 to 




7r/90 


tt/5 


0.45 














100 exp 


-t] 




SLOW.START.3'^ 


1.0 


5.0 


0.0 


0.75 


0.01 


20 to 
100 exp 


-t] 


7r/90 


tt/5 


0.55 



Note. — A representative sample of our non-systematic exploration of parameter space; not 
all runs attempted are included. "Interpolation" runs refer to those runs which are a natural 
extension of Tables ^ and [2] whereas "other" runs include substantially different atmospheres and 
jet behaviors. Sec text in Section 3.2 for discussion. 



^ This is the standard run with axisymmetric Kolmogorov spectrum density perturbations intro- 
duced to the atmosphere. 

^ Shells of material with a peak amplitude p — 3.0 were superimposed on the standard atmosphere; 
shells are generated along the 2-axis by using 2D sine or sine functions with a period of 2 ~ 0.5. 
Djct = 2OC3 from t = 0.0 to t = 0.25 and ujot = 100exp[-(i - 0.25)] thereafter. 

neither case are wings promoted, as each successive brief 
outburst deposits most of its thrust at the ends of the 
lobes, far from the midplane. The chief reason why such 
intermittency does not produce prominent wings is that 
wings expand subsonically. Since the duty cycles of the 
intermittent jets are much shorter than the crossing time, 
the cocoon expands as if a moderately powerful jet of 
constant velocity were powering it. Unlike in single de- 
caying outbursts, the intermittent jets deposit most of 
their thrust far from the midplane. On the other hand, 
intermittency on the timescale of an e-folding time of a 
decaying jet can be effective. 

4. SIMULATION LIMITATIONS 

Because we have only solved the equations of hydro- 
dynamics in evolving our simple models, it is worth con- 
sidering the impact of additional complexity. 



4.1. Missing Physics 

Our models do not include magnetic fields, special rel- 
ativity, or radiative losses, and our simple setup does not 
take into account complex jet or atmospheric structure, 
feedback, and other processes that may be important to 
radio galaxy morphology. In practice, we do not believe 



these omissions invalidate our results. For example, spe- 
cial relativity and magnetic fields must play crucial roles 
in determining the character and transverse structure of 
the jet, but we only insist that our simulated jets repro- 
duce the coUimation, hot spots, and back-fiows of FR II 
sources. In other respects (e.g. radiative efficiency), the 
jet and AGN are inside a "black box." Because we are 
concerned with the behavior of lobe material, we believe 
these omissions are justified. 

Likewise, radiative losses (depleting of lobe energy via 
synchrotron emission and inverse Compton scattering) 
only become important relative to adiabatic losses at late 
times and are not in a position to infiuence the formation 
of wings. Radiative losses from the ICM are also irrel- 
evant on the timescales of AGN outbursts, so feedback 
(the connection of jet power to the amount of material 
crossing the inner boundary) resulting from cluster cool- 
ing is unimpor tant on these timescales. As to feedback 
from backflow ijAntonuccio-Delogu fc SiiSl2010l ). there is 
no obvious recipe to describe what happens to the mass 
flowing across the inner boundary, since the inner bound- 
ary radius is much larger than the nuclear engine, but as 
the back-flows straighten, less material will have veloc- 
ity vectors pointing towards the origin and eventually 
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the AGN would cease to be "fed" by backflow. This is 
consistent with our entirely artificial recipe for decaying 
jets. 

On the other hand, magnetic fields (that must be in- 
jected with the jet) and relativistic jets (which have 
a higher thrust for a given mass-energy density) may 
strongly influence the behavior of lobe material. 

Toroidal magnetic fields may help collimate the back - 
flow and retard its mixing (e.g. iBraithwaitd |2Q10| ). 
The apparent continuity of fields in some wings (e.g. 
NGC 326 in Murgia et al. 2001) is suggestive, especially 
considering that in our models the lack of coUimation 
leads to wide wings. However, iHuarte-Espinosa et al.l 
(|2010D find in simulations of FR II sources that turbu- 
lence sets in within the cocoon even when ordered fields 
existed earlier, so it is unclear that the magnetic fields 
in the wings are actually toroidal. As to the realism of 
the jet, it is not our goal to understand the jet physics in 
detail, but we note the persistence of experimental hy - 
personic fluid jets in recent work bv lBelan et al.l ()2011l ). 
which the authors use to argue that magnetic collimation 
is only required near the base of the jet. 

The importance of r elativistic jets to t he cocooon 
morphology is unclear. IKomissarov &: Falld (1996) find 
that, for jets matched by velocity, pressure, radius and 
power, relativistic jets produce preferentially wider co- 
coons than classical ones. The authors suggest that 
this effect may be accounted for b y the higher thrust 
in the relativistic case. However, IRosen et al.l ()1999f ) 
find the opposite when using the same matching condi- 
tions: relativistic jets tend to produce narrower cocoons 
compared to classical jets (more specifically, jets with 
a higher Lorentz factor have smaller cocoons). With- 
out matched relativistic simulations it is impossible to 
directly test the importance of speci al relativity on our 
coc oon sizes, but it is reinarkab le that lRosen et al.l (|1999[ ) 
and lKomissarov fc Fall^ (|1996D agree that the qualitative 
nature of the cocoon is unchanged. Including relativity 
may modify W/L but would not fundamentally change 
the appearance of our winged sources because relativistic 
physics in the lobes would not change the processes which 
produce wings. In other words, we believe that special 
relativity is a second-order effect akin to magnetic fields. 

4.2. Initial Conditions 

It is worth asking whether wings form only be- 
cause our atmospheres are "smooth" , i.e. because we 
have not included turbulence in the ICM or struc- 
ture near the center (where the environment is pre- 
sumably complex). This question is particularly impor- 
tant given the obser ved int eractions between jets and 
molecular clouds fe.g. ILv et al. 2005), stellar shells (e.g. 
iGopal- Krishna fc Wiit all 20101 . gas in companion galaxies 
(jEvans et al. 2008.) . and more generally th e complex nu- 
clear environments of active galaxies (e.g. iRosario et al.l 
|2010[) . To this end, we have introduced density pertur- 
bations mimicking stellar shells along the major axis at 
various intervals, tested the impact of a single ring near 
the nucleus, and introduced Kolmogorov-spectrum tur- 
bulence to the ICM. None of these structures alter the 
same basic evolution of winged sources. Additional jet 
physics would need to be in place in order to determine 
whether the jet itself can be b ent by in teraction with 
high density pockets (|GoDal-Krishna fc Wiitai2010i ). All 



of these tests use density perturbations with a maximum 
amplitude of the core density pQ. Finally, bulk flows in 
clust ers are clearly im portant to radio galaxy structure 
(Mor sonv et al.|[2010[ ). although most XRGs are large, 
strongly bridged sources in which this may not be a defin- 
ing effect. 

4.3. Artifacts 

It is encouraging that our simulated radio galaxies re- 
produce the basic features of other simulations in the 
literature (e.g. Reynolds et al. 2002 ; Heinz et al. 200^ 
Vernaleo fc Reynolds! 120071 : lAntonuccio-Delogu fc SilH 
20101) . but it is worth considering the impact of bound- 
ary conditions and the jet structure on the morphology 
of the lobes — our setup is designed to produce lobes like 
those of radio galaxies with physics and conditions in- 
herently different from those encountered in nature. The 
mixing experienced by the lobes is also artificial (natu- 
rally set by the size of the grid zones, which vary along r 
and 6), but mixing is too slow a process to suppress the 
formation of wings. 

The inner boundary is the most significant artifact of 
our simulation because it represents no physical analog 
but rather allows us to hide the AGN and the jet col- 
limation mechanism. In addition, backflow crossing the 
boundary slows down and disappears. This has two con- 
sequences: not all the backflow can be harnessed, and 
flows may be effectively directed around the boundary 
by eddies as slowing backflow crosses the boundary. In 
the flrst case, the amount of material is to small to in- 
fluence the wings. In the second, we flnd that the inner 
boundary must be small compared to rg (keeping the 
physical size of the jets flxed) in order to prevent the ini- 
tial cocoon from being unduly influenced by the sphere. 
Inner boundary spheres which are too large tend to pro- 
mote material flowing around them, hence promoting 
early wings, but suppress midplane mergers and hence 
the prominence of later wings. The initial interstices 
also survive for longer for oversized Tinner ■ For typical 
atmospheres (0.5 < rg < 1.0), rjnner must be somewhat 
smaller than 0.1 code units; Hnncr — 0.05 is a compromise 
between a small impact on the initial cocoon and reason- 
ably sized grid zones. Runs with Hnner = 0.01 do not ap- 
pear substantially different from those with Hnncr — 0.05, 
and the timestep is unreasonably small. 

We precess the jet injection footpoints rapidly (207r Hz) 
in a very tight circle around the poles in order to break up 
the jet symmetry and spread its thrust out over a larger 
working surface. We only require that our jet reproduce 
internal features seen in other hydrodynamic simulations 
along with the terminal shocks that give rise to the back- 
flows and lobes. Without breaking up the jets, the jet 
hea d travels very quickl y and p roduces unrealistic lobes 
(e.g. lVernaleo fc Revnoldi[2006[ ): in MHD simulations we 
would expect helical instabilities to break up the jet. Our 
precession scheme is thus an artifice that is motivated by 
observables but whose recipe is not an at tempt to model 
a physical process fc.f. IHeinz et al.l [ 20 06t). However, one 
might worry about its influence on backflow and wings 
because the precession angle is a free parameter (subject 
to the condition that a < Ojet)- In other words, we can 
tune the width of the lobes and the rate of growth of the 
radio galaxy within some narrow range of parameters. 
Since stalling the jet head contributes to wing promi- 
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nence via shortening the active lobes, is the axial ratio 
W/L artificially high? We do not believe so. Even with 
tiny precession angles, we form long wings in highly ec- 
centric atmospheres (the only atmospheres where long 
wings form). Moreover, the jet head advance speed also 
depends on the kinetic luminosity as a function of time. 

5. PROPERTIES AND PREDICTIONS OF THE BACKFLOW 

MODEL 

We have successfully produced winged and X-shaped 
sources solely by redirecting back-flow with static en- 
vironmental pressure gradients, and are now in a 
position to compare our simulated sources to ob- 
served XRGs (Figure [TT|) in order to critique pro- 
posed wing formation mechanisms involving back-flows. 
In these models, wings are produced as plasma flow- 
ing back from the jet heads is deflected into a di- 
rection misaligned with the jets. How thi s de- 
flection occur s is u n certain ([Leahv fc "W illiams 119841: 
Worrah e t all Il995t iCapetti et alj l2002u .Kraft et al.l 
2005; Go pal-Krishna fc Wiitall2010D . and until now the 
backflow model has not been rigorously investigated with 
three-dimensional simulations. 

Here we outline a series of predictions (expectations) 
for properties which are, at least in principle, observ- 
able if our simulations are an adequate representation 
of wing formation. These predictions are also summa- 
rized in Table We then assess the backflow model in 
light of these predictions and briefly discuss the implica- 
tions of our work for the wider sample of distortions to 
the canonical double-lobed morphology of bridged radio 
galaxies. 

5.1. Predictions 

X-shaped and wmged sources are one family — The most 
natural consequence of the backflow model is that short 
and long wings are produced in the same way. Short 
wings should then be more common than long ones since 
the long wings are harder to produce. Although this is 
an obvious point, we note it as a potential observational 
test because sources with shorter wing s and the can- 
didate X-shaped sources (jCheungl 120071 ) have not been 
tested for several of the trends seen in X RGs (most im- 
portantly, the jet-major axis correlat ion ICapetti et al.l 
[2002: Sa ripalli &: SubrahmanvanI |2009() . Our simulated 
sources predict that winged sources will fall in line, al- 
beit with a greater spread in parameter space. 

Projection almost always enhances W/L — Wings in our 
models (even those produced by triaxial atmospheres) 
tend to be wider than the lobes and hold their shape 
better in rotation and projection. Hence, we predict that 
projection almost always enhances W/L. In other words, 
some observed sources with high VF/Lobs probably have 
an intrinsically smaller aspect ratio (e.g. Figure [TT|) . If 
the projection angle can be worked out for a sufhcient 
number of XRGs and winged sources, we expect the lobes 
in a number of XRGs to be shortened via projection, 
i.e. the intrinsic distribution of W/L is shifted from the 
observed distribution (an idea of the observed distribu - 
tion may be found in lSaripalli fc SubrahmanvanI I2OO90 . 
If this does not turn out to be the case, the backflow 
model would need a collimation mechanism (see subsonic 
expansion below) to be consistent with observations. Un- 



fortunately, this prediction makes it difficult to quantita- 
tively compare the simulated population to the observed 
on because in most XRGs the projection angle is un- 
known. 

Long wings require high ellipticity — Our simulations pre- 
dict that intrinsically long wings can only be produced 
in very elliptical atmospheres (Figure [7]). The lowest e 
for which we can make convincing wings is e ~ 0.45, 
and very long wings require e > 0.55. These values are 
much higher t han the average powerful radio galaxy host 
(epcak ^ 0.2; iSmith fc HeckmanI Il989f ) and even most 
XRGs. A high initial atmospheric ellipticity (or higher 
order asymmetry) is a clear prediction of our models, but 
confirming this behavior in observed sources is difficult. 

For instance, a few sources with very small host e and 
large W/L (the best example is NGC 326) appear to 
have wings originating outside the ISM, so it is not ev- 
ident that the ellipticity of the ISM is always the rele- 
vant value for comparison. Further, the ellipticity of the 
ISM may differ su bstantially from t hat o f the host galaxy 
on small scales (jPiehl fc Statled 12007^ eve n when in 
broad agreement on galactic scales (jHodgc s-Kluck et all 
[2OT0h . In this case, the best example in the literature is 
3C 403 (Figure [H where eisM = 0.059 and eopticai = 0.25 
(|Kraft et al.ll2005D . Indeed, we find that elliptical atmo- 
spheres of moderate e with embedded disks are as effec- 
tive at producing wings as single-component atmospheres 
with extreme e, but the disks are largely destroyed by the 
radio galaxy. High resolution imaging may be required 
to see any disk remnants. The triaxiality of elliptical 
galaxies also plays a role, since we may not measure a 
true ellipticity. Finally, because projection tends to en- 
hance W/L, extreme values of e 0.75 (as we use in 
the standard atmosphere) are not required to produce 
observed long wings. All else being equal, we still ex- 
pect the intrinsic W/L to be correlated with e, but the 
caveats outlined above make it a difficult proposal to test 
presently. 

Wings require a jet pointed nearly along the major axis 
of an anisotropic environment {A PA ~ Oj — As ex- 
pected (based on the studies in ICapetti et al.l 120021: 
Saripalli fc SubrahmanvanI 120091 : iHodges-Kluck et all 
20101) . wings require a strongly asymmetric environment 
in which the jet is stalled by progression along the long 
axis (Figure ^ . The wings grow in the favorable pres- 
sure gradients along the minor axis or axes. Our models 
expect a fairly strict requirement for the jet to be within 
~ 15° of the major axis for substantial wings to be pro- 
duced via the overpressured cocoon channel (the usual 
method for seeding wings in our models). The require- 
ment is even more stringent in the event that a disk is 
present, indicating that disks are likely to be important 
only in a minority of cases. 

As APA increases, the wings become shorter (FigurejT]) 
as well as increasingly associated with the lobe on the 
same side of the major axis. These lobe-wing pairs thus 
have mirror symmetry about the midplane. This is in 
agreement with radio maps of such sources. Since wings 
long enough to qualify the source as "winged" can be pro- 
duced for a 15° range of major axis-jet separation and 
truly X-shaped sources are only (intrinsically) produced 
for APA < 5°, we would expect ~ 2/3 of winged sources 
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Table 4 

Model Predictions 





Expectation 


Observed? 


Reference 


1 


OliUi L W illga IjL'llilllUii, t iilLl illolLclll V J ItJlig 


iNtj 111 oouy odiiijjit 


SS09 






I iTi/^lii/HiTirr i M/ri 




z. 


Projection tends to enhance W/-L 


Undetermined 




Q 

o. 


W/L correlated with cism 


Sample size too small 


UTT^ 1 n 

rllviU 




(weaker correlation with host e) 


INo : 




4 




Yes 


Pn9 SSDQ HKl n 




Llit: llidJUi dJtiO 






r: 
0. 


enhanced by higher pressure, 




Tin 




small atmospheres 


Vat3? 

leS! 


oouy 


6. 


Most XRGs should be weak FR lis 


Yes 


COQ Lin 




(and hence strongly bridged sources) 


Yes 


LW84 


7. 


Intermittency on scales of 3-10 Myr 


Undetermined 






produces longer wings 






8. 


Wings are fainter than lobes 


Yes 


LW84 


9. 


Wings have flat spectral indices 


Mixed results 


LR07 


10. 


Backflow follows existing channels 


Yes? 


SS09 




(collimation requires existing channels) 






11. 


Wings grow subsonically 


No? (large wings) 






(hence wing AGN outbursts are old) 


No? 


DT02 


12. 


Flow speeds in wings should be 


Undetermined 






transonic for lobe material 








Backflow model piggybacks on other hydrodynamic 


Yes? 


ssoa 



models 



References. 



— C02: 



ICapetti et al.l ll2002f): C07: ICheund (1200711 : 09: ICheung et al.l ll2009f) : 

DT02: | Dennett-Thorpe et al'l I I2002I): HKIO : IHodges-Kluck et al.l ||2010I)Tl10: [ Landt et al.l 1 12010) : LW84: 
ILeahv fc Williamj ngS^) : LR07: ILal fe Raol 120071) : SS09: ISaripalli fc SubrahmanvanI II2009I) 



Note. — Predictions for XRG properties if our model is accurate; note that these may not hold for the 
"backflow model" generally as phrased in prior work. Predictions correspond to arguments made in Section 4. 



to have lobe- wing pairs, with acute angles between a lobe 
and its wing produced for smaller angles and obtuse an- 
gles for larger ones. Thus, we suggest that the existence 
of a large number of sources with lobe- wing pairs rather 
than true X-shaped morphology is consistent with the 
backflow model. However, we would also expect that 
wing length would decrease with increasing APA; this 
has not yet been measured. 

Longer wings are produced in higher pressure, smaller at- 
mospheres — Steeper pressure gradients occur for smaller 
core radius and higher core density/pressure (Equa- 
tion 12). Since the steepness of the pressure gradient 
influences the rate at which wings grow, these small, 
dense atmospheres (e.g. the ISM) are better at produc- 
ing wings. Moreover, wings can only form via the over- 
pressured cocoon channel if the cocoon can escape the 
central core region and its high density before it comes 
into pressure equilibrium. The top panels of Figure 12.21 
demonstrate this point (see also Figure [7]). Finally, a 
high core density also resists the jet advance along the 
major axis, allowing the wings to grow longer. Thus, our 
models expect that the hosts of XRGs (where the jet is 
presumably pointed along the major axis) have higher 
pressure, on average, than those no rmal radio galaxies 
with jet geometry favorable to wings. iLandt et al.l (poTof ) 
find that the nuclear regions of XRGs have high temper- 
atures (T ^ 15000 K) indicating that these regions may 
be overpressured. Although this does not directly corre- 
spond to our models of relaxed, isothermal atmospheres, 
our requirement for a small, high pressure environment 
is in agreement with their work. 

Most XRGs should be weak FR lis — Another key ingredi- 
ent in our models is a decaying jet that begins as a power- 
ful FR II source and decays to luminosities more typical 



of FR I sources. A powerful jet is required to produce the 
overpressured cocoon and drive wings early on, whereas 
once wings begin expanding subsonically, a weakening 
jet (which advances increasingly slowly along the major 
axis of the atmosphere) allows the wings to become quite 
prominent. If the atmosphere is much smaller than the 
radio galaxy (as in many cases) , a decaying jet will still 
grow more and more slowly, allowing subsonically grow- 
ing wings to keep up. This finding naturally explains the 
observations that XRG radio powers tend to lie near the 
FR I/II "break" (jCheung et al.ll2009f ) while also possess- 
ing strong bridges associat ed with powerful FR II sources 
(Leahy & Williams 1984|). ILandt et al.) (I2010D argue that 
XRGs are the archetypal transition population between 
FR I and FR II sources, with about half the XRGs in 
their sample having weak emission lines from the AGN 
(weak-lined FR lis are otherwise uncommon) . This sug- 
gests that there is indeed an evolutionary progression 
in the AGN and jet luminosity of XRGs. Our models 
predict a fast-rise exponential-decay profile in which the 
XRG spends most of its lifetime a s a weak FR II source 
(Figure [8]). However, we note that iBesd (|2009t ) find in an 
SDSS study tha t FR I and II sources are not as obviously 
separated as in iLedlow fc OwctI ([19961 ) and that FR II 
morphology occurs for a variety of radio powers and host 
galaxy masses. 

Intermittency is another important prediction of our 
models for observed XRGs, but the action of intermit- 
tency in our simulations may represent a more complex 
underlying process. Intermittency on timescales similar 
to the e-folding time of a decaying jet (near the early 
passive phase) is effective in our simulations because it 
allows the jet to bypass the overpressured cocoon stage 
and form a bow shock within the radio lobe itself, thereby 
depositing a large amount of backflow into the wings from 
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a very powerful jet close to the wing bases. Reinjection is 
effective in a relatively narrow window of time (3-10 Myr 
depending on the size of the radio galaxy, e.g. Figure [5]) 
between the active and passive phases, bu t this may be 
plausible: the ripples in the Perseus cl uster (Fabia n et al.l 
[2006) and Abell 2052 (Bla nton et al.l , 2009) do imply an 
AGN duty cycle of about 10 Myr. This is consistent 
with our simulated sources between 50-100 kpc across. 
Thus, it is conceivable that regular intermittency is a vi- 
able mechanism for enhancing wings, although we note 
that if the intermittency is very rapid or very slow, it 
is ineffective (Section 3.2). On the other hand, intermit- 
tency could be replaced in our models by any mechanism 
which drives back- flows nearer the base of the wings such 
as the motion of denser ISM into the path of the jet. We 
note that intermittency may also be effective for wings 
attached to the primary lobes farther from the nucleus; a 
reinjected jet quickly re-establishes the cocoon and drives 
strong back-flows along its length as it rapidly reaches 
the prior hot spot. If the jet reinjected is similar to the 
outburst which initially formed the wings, we would still 
expect XRGs mostly to be weak FR II sources. 

Wings should be fainter than the primary lobes — The ac- 
tive lobes of XRGs are typically brighter than the wings. 
This behavior persists at lower (MHz) frequencies and 
is therefore not obviously attributable to spectral age- 
ing, alt hough spectral a geing may be in play for some 
XRGs (|Lal fc Radl2007| ). Our models expect the wings 
to be fainter because they are substantially wider than 
the collimated primary lobes yet contain (in most cases) 
less lobe material. We therefore suppose that the wing 
material spreads out and decreases the magnetic field en- 
ergy density Ub- 

Assuming that Ub decreases with increasing volume 
(i.e. assuming something like equipartition conditions, 
although precise equipartition is not required), the syn- 
chrotron emissivity in the wings will be substantially 
smaller than that in the lobes. In other words, "spheri- 
cal" wings a factor of ^ 2 wider than the active lobes will 
have a factor of 8 smaller Ub than in the lobes. In op- 
tically thin conditions (a good assumption), the surface 
brightness obtained from integrating through the wings 
only recovers a factor of 2 in the wings, so we expect 
them to be a factor of ^ 4 dimmer than the active lobes 
(more generally, the square of the ratio between the wing 
width and the lobe width) for electrons of the same 7. Of 
course, in projection the wings may appear even dimmer 
since the lobes will tend to be shortened and therefore 
increase in surface brightness (ignoring relativistic dim- 
ming of the counterjet) whereas the wings will not. 

Clearly, this prediction rests on a number of assump- 
tions not included in our models and is therefore some- 
what weak. Since the wings are filled with turbulent 
plasma, it is also possible that the magnetic field strength 
is increased by winding up of fields as the material mixes. 
Moreover, the jet power (and therefore the back-flow 
speed) is time dependent in our simulations, so the ra- 
tio between the surface brightness of the lobes and the 
wings would be time dependent in our model as well. 
Nonetheless, the diffusion of material in the wings and 
the aforementioned projection effects seem likely to dim 
the wings in the backflow model. 



Backflow follows existing channels — Most of the wings 
in our simulated sources are seeded by the pressure- 
driven expansion of an overpressured cocoon early in the 
source's life. This cocoon expands asymmetrically be- 
cause of the asymmetric pressure gradients and forms 
proto-wings which are later bolstered by the merger of 
laminar back-flows near the midplane sending material 
into the wings during the active phase. The laminar 
back-flows themselves are ineffective at drilling new chan- 
nels. During the overpressured phase they acquire a vor- 
ticity near the terminal shocks that leads them to follow 
the contact discontinuity of the cocoon and ultimately 
flow back towards the AGN (Figure [3]). During the ac- 
tive phase, they are mostly straight and flow towards the 
midplane. 

Lobe material is also ineffective at producing channels 
because it follows the path of least resistance: in a con- 
flned ellipsoidal cocoon, it will simply spread out to in- 
crease the pressure throughout the cocoon rather than 
break out in a particular direction. This can be seen in 
the false synchrotron maps in Figure 12.21 These maps 
assume equipartition and are not indicative of what the 
sources would look like in the GHz bands. Rather, they 
effectively trace pressure in the cocoon, and it is easy to 
see that the pressure in the confined cocoon on the left 
side of the figure remains high relative to the X-shaped 
source (the false radio maps use the same scale). 

Apart from proto-wings produced by the early 
pressure-driven expansions, Kelvin- Helmholtz and 
Rayleigh- Taylor instabilities produce channels which 
backflow reinforces. Given sufficient time and a powerful 
jet, these whorls and fingers can become wings in their 
own right, but likewise this is due to the growth of 
the instabilities and buoyancy rather than redirected 
laminar back-flows. 

Notably, some sources exhibit wings which do not em- 
anate from the center or are unlikely to have been pro- 
duced via an overpressured cocoon (e.g. NGC 326). Be- 
cause it is difficult to channel the back-flows, we hypoth- 
esize that if the backflow model is correct, some seed 
proto-wing was necessary to produce such sources. The 
origin of these proto-wings is unclear. 

Wings expand subsonically — Related to the expectation 
that back-flows are not directly responsible for drilling 
channels is the prediction that the wings of radio galaxies 
expand subsonically for most of their lifeti mes ( Figures [5] 
andffl) . This expectation is contrar y to the | Cap etti et al] 
(I2OO2D "overpressured cocoo n" and lWorrallet al., (|1995| [ 
or Le ahy fc Williams! ()1984[ ) "buoyant backflow" propos- 
als in the literature (Section 1) and poses a serious prob- 
lem for the backflow model in terms of explaining long 
wings. We therefore discuss this point in some detail 
below. 

The problem of subsonic expansion can be broken into 
two distinct objections: (1) If the wings grow subsoni- 
cally while the primary lobes grow supersonically, how 
can wings be longer than the lobes? and (2) Subsonic 
expansion implies wing lifetimes of more than 100 Myr 
for sources hundreds of kpc across. The first objection is 
easier to reconcile with our models because our jets pro- 
vide a natural mechanism for the subsonic expansion of 
the lobes, and it is not clear whether most radio galaxies 
(even powerful ones) really do expand very supersoni- 
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Figure 9. Gallery of false equipartition synchrotron maps (j^ oc pj p'^'^) showing the resemblance between our winged sources 
and observed XRGs, and the importance of projection. Left: Face-on view showing intrinsic structure. The red ellipses show 
the core radius isobar of the model atmosphere. Center: Rotated and projected view. Right: XRG analog (various resolutions). 
We only show a few examples of simulations and note that any one simulation can reproduce several observed sources; more 
than 60 of the sources in Cheung, (2007) can be plausibly reproduced by our simulations (in terms of appearance). 



cally. Further, as we have noted, projection tends to 
enhance wings. However, the lifetimes of large sources 
are more difficult to explain. 

To see this, consider that for a rich cluster atmo- 
sphere with kT ^ 10 keV (this corresponds to a '--^ 
1000 km s~^ or, conveniently, 1.0 kpc Myr~^), subsonic 
growth implies that wings 100 kpc from base-to-tip would 
be at least 100 Myr old. In group environments, where 
Cs < 500 km s~^ an d monotonically declines at large 
radii (jSun et al.ll2009D . the longest observed wings (e.g. 
3C 315 whose wings span ^ 400 kpc) could be almost 
500 Myr old! Even if AGN activity persisted this long, 
synchrotron cooling might set in: if the bulk flows of re- 
plenishing backflow travel at a lobe sound speed lOc^, 
they would take between 20-100 Myr to replenish the 
leading edges of the wings. Depending on the break fre- 
quency and assuming magnetic fields in equipartition {B 
of a few jiG) and a Lorentz factor 7 of severalxlO^, this 
may easily exceed the radiative cooling time. 

To solve this problem, variations of the backfiow model 
require that the wings actually do expand supersonically. 
In the "overpressured cocoon" model ()Capetti et al.l 



[200l . supersonic wing growth is achieved through sus- 
tained pressure-driven expansion, whereas i n the "buoy- 
ant backflow" model (jWorrall et al.l Il995j ). wmgs are 
driven by the nearly free expansion of diverted hyper- 
sonic back-flows. Our models reproduce overpressured 
cocoons and hypersonic back-flows, but the overall wing 
advance speed is subsonic. 

For the overpressured cocoon model to produce long 
wings, a sustained overpressured state must be main- 
tained to drive outfiows. However, in our models the 
highly overpressured cocoon that forms soon after ig- 
nition quickly expands to reach pressure equilibrium: 
the supersonic expansion phase lasts at most around 
15% of the lifetime of the radio galaxy regardless of 
the jet we inject (e.g. Figure |4]). As the radius of the 
(isobaric) cocoon grows, its pressure falls much more 
rapidly than does that of t he surrounding at mosphere 
at a similar radius (see also iZanni "eFan[2003L Fi gure Z 
in their paper), so the advance speed of the lateral co- 
coon expansion falls precipitously from a peak of about 
Mach 2. Thus, overpressured expansion appears to be 
relevant only very early in the life of the source and by 
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Figure 10. Figure 9 continued. 



itself can produce only short wings. Indee d, this con- 
clusio n is supported by earlier simulations (jZanni et al.l 
I2003f ) where even in advantageous triaxial atmospheres, 
the overpressured cocoon produced at most an intrinsic 
W/L of ^ 0.5 (the jets in their work were held at con- 
stant velocity and W/L would th us decrease over time). 
iSaripalli fc SubrahmanyanI (|2009[ ) suggest that the over- 
pressured state may instead result from backflow piling 
up upon reaching the central (dense) region of the galaxy. 
Our models do not support this scenario, since the co- 
coon is contiguous and nearly isobaric throughout the 
active phase. 

The buoyant backflow model argues that a combi- 
nation of buoyancy forces and wings driven by hyper- 
sonic flows produc e large X- s haped sources. For in- 
stance, Dennett-T horpe et al.l (|2002() estimate that the 
current outburst in 3C 403 (Figure [1]) started 16 Myr 
ago. If this outburst was solely responsible for gener- 
ating th e ^ 100 kpc wing s hydrodynamically (as pre- 
ferred bv lKraft et atll2005f) . the average expansion speed 
of the wings would have to be ^ 8000 km s~^ (0.027c). 
In contrast, a typical sound speed for galaxy groups is 
Cs 5, 500 km s"^, a factor of 16 smaller. Likewise, 
iWorrall et al.l (jl995D prefer hypersonic expansion. What 
could drive such wings? The answer invoked by these 
authors is redirected back-flows, which can be accel- 



erated at the terminal hot spots up to a few percent 
of c (based on the spect ral ageing-distance method of 
[Alex ander &: Leahvlll987l ). The observational inference 
of high bac kflow speeds is consistent with our models and 
others fe.g. lAntonuccio-Delogu fc Silk|[2010t ). which find 
back-flow speeds of about t wice the lobe sound speed 
Cs lobe = lOcs). Further, ISaripalli fc SubrahmanvanI 
2001 argue that the collimated morphologies of some 
wings requires some driving force. 

However, it is not evident that these laminar back-flows 
can be harnessed to drive wing expansion (see the prior 
prediction). Subsonic wing growth in our simulations is 
a natural consequence of the tendency of backflow (lobe 
material) to mix and expand. The fast back- flows from 
the terminal shocks merge near the midplane, dissipat- 
ing their speed and driving flows into the wings which are 
subsonic relative to the low-density material in the lobes 
(Figure [3]) , although they are still supersonic relative to 
the background. These flows then expand and decelerate 
further in the wings, dissipating their thrust over a very 
large area. Hence, detection of 10'^ — 10* km s~^ flows 
in the wings would be insufficient evidence for hyper- 
sonic wing expansion. Supersonic wing expansion would 
require coUimation to prevent thrust dissipation. 

Can a coUimation mechanism be found? An obvi- 
ous possibility we have not included is ordered magnetic 
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fields, since particles c an move along them much more 
easily than across, but (jHuarte-Espinosa et al.l[2010l ) find 
that in FR II sources, lobe magnetic fields become tur- 
bulent on timescales of ~ 10 Myr. More generally, we 
note that (weak) shocks found around double-lobed ra- 
dio sources often imply cocoon expansion near Mach 1. If 
a powerful jet cannot drive highly supersonic expansion, 
it seems unlikely that less powerful, more disorganized 
backflow could do so. 

Lastly, it is important to emphasize that subsonic ex- 
pansion is only a serious problem for very large wings or 
very young AGN with a large W/L ratio. Our simula- 
tions can reproduce most winged and X-shaped sources 
with subsonic expansion if the AGN outburst is between 
10-100 Myr old instead of 1-10 Myr. These simulations 
are also not the final word, as they only solve the equa- 
tions of hydrodynamics in a relaxed, non-dynamic atmo- 
sphere. 

5.2. Assessment 

Our results, taken in context, strongly implicate a hy- 
drodynamic origin for X-shaped sources and a common 
origin for winged and X-shaped sources. We can produce 
bona fide XRGs (i.e. with intrinsic W/L > 0.8) with a 
single outburst from a jet with a plausible time depen- 
dent kinetic luminosity on timescales broadly consistent 



with observed sources. The radio galaxies we produce 
are consistent with prior simulations, and faithful repro- 
ductions of observed radio galaxies can be generated by 
tuning the viewing angle (Figure [TT|) . These sources are 
also consistent with most observed properties of the XRG 
population outlined above, and our simulations repro- 
duce elements of the overpressured cocoon and buoyant 
backfiow models while relying solely on the interaction of 
the radio lobe material with an anisotropic environment. 
Hence, the backfiow model remains a strong contender 
for the origin of XRGs — it would be surprising if some 
XRGs were not produced in this manner. 

However, our simulations also present new challenges 
to the proposals in the literature: we find that long 
wings require unusually elliptical atmospheres and ex- 
pand subsonically, making very long wings difficult to 
explain. Having examined the deficiencies of our simu- 
lations, we cannot identify an obvious internal remedy; 
we note that pre-existing channels are necessary (in the 
backfiow model) for some sources in the literature and 
may be generally important. 

The final prediction of our simulations is therefore that 
the XRG population is heterogenous and the backflow 
model works in tandem with other hydrodynamic models 
rather than solely on its own. In short, this is because 
our models make it too difficult to produce XRGs given 
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their observed frequency. 

Consider that winged and X-shaped sources make up 
5-10% of double-lobed radio sources. If, as in our mod- 
els, X-shaped sources are produced by a fortunate coinci- 
dence of jet geometry, ISM cllipticity, and jet power, we 
might imagine that the fraction of X-shaped sources is 
some function of each of these variables. In the simplest 
case, where each of these factors contribute indepen- 
dently to promoting wings, the fraction of XRGs might 
look like 

/XRG CX f^Af^fL.L^ (14) 

where /xrg is the fraction of double-lobed sources which 
are X-shaped and the other / values represent the frac- 
tion of sources for each variable which meet the threshold 
criterion for wings. The exponents are unknown mea- 
sures of the relative importance of each variable. Ob- 
viously, the fraction of XRGs must be smaller than the 
fraction of sources which meet any one criterion. 

Now, consider that virtually all double-lobed radio 
galaxies emanate from elliptical hosts, and have appar- 
ently random jet-major axis orientations (as opposed 
to the very weakest radio-emitting AGN whose jets 
seem to be pointed alon g the minor axis of their hosts; 
iBrowne fc Battvd |2010[ ). Our models suggest that a 
jet within ~ 15° of the major axis of its host is re- 
quired to produce wings. In a uniform distribution of 
radio jet position angles, we would thus expect that 
/pa '~ 15/90 PS 0.17. Assuming all powerful jets have a 
kinetic luminosity function conducive to forming XRGs, 
we still require about half of all sources where the jet 
is co-aligned with the major axis to have high ellipticity 
(assuming the exponents are all equal to unity). Given 
that ellipticity is clearly a very sensitive parameter and 
that the peak e of po werful radio galaxy hosts is far below 
the wing threshold (|Smith fc Heckmanlll989| ). it is clear 
that our simulations (in this oversimplified formulation) 
underpredict the observed frequency of winged sources. 
In other words, wings are too difficult to make in our 
models. 

To reconcile this result with the observations, either 
our models must be fundamentally deficient with respect 
to the behavior of backflow or they require alternate 
mechanisms to form proto-wing channels such as the jet- 
stellar shell interaction model of iGopal-Krishna fc Wiital 
(|2010( ) or the jet-merging ISM explanation for Z-shaped 
sources in Zier (2005). As we have seen, backflowing ma- 
terial reinforces any proto-wings and can turn them into 
full-fledged wings if the original mechanism fails to do 
so; making the channels in the first place is what the 
backflow model cannot easily do. 

Thus, we propose that the backflow model has a com- 
mensal relationship with other hydrodynamic mecha- 
nisms for forming wings by reinforcing and growing any 
pre-existing channels accessible to a jet pointed along 
the major axis of its host galaxy. These channels would 
naturally grow most easily along the steepest pressure 
gradient. In cases where the jet power and atmosphere 
size match appropriately, these channels would be pro- 
duced by the expansion of an overpressured cocoon as in 
our simulations. In other cases, channels could be pro- 
duced by the interaction of the jet with structure in the 
ISM or by minor mergers. In this scenario, we would still 
expect to see most of the predictions of our simulations 



listed above for the case of pure backflow in a relaxed 
atmosphere. 

5.3. Distortions to the Canonical Double-Lobed FR II 
Radio Galaxy 

At this point, it is worth r evisiting the buoya nt back- 
flow model as phrased in iLeahv fc WilliarasI (p.984), 
where XRGs are unifled with other distortions to the 
canonical FR II model (about two thirds of strongly 
bridged radio galaxies show central distortions) by 
the deflection of backflow around a denser medium. 
These distortions may also be related to radio galax- 
ies with interrupted bridges b y the "superdisk" model 
(jGopal-Krishna fc Wiital |2009() . in which gas displaced 
by a galaxy merger is believed responsible for docking 
the lobes in an asymmetric manner. 

Because the backflow model as expressed in our simula- 
tions depends only on a fortuitous combination of atmo- 
sphere morphology and jet behavior to produce winged 
sources, we can indeed reproduce the basic bridge distor- 
tions of lLeahv~fc Williams (1984 ). However, we cannot 
reproduce long Z-shaped sources and some asymmetric 
distortions which do not appear to be the result of bulk 
motion or cluster turbulence. These include sources with 
a single wing on one side, sources with wings which them- 
selves bend dramatically, sour ces with one FR I lobe 
and one FR II lobe (HYMORs; iGopal-Krishna fc Wiital 
I2OOOI) . A gain, the main hindrance seems to be the inabil- 
ity to channel backflow; jets which bend may solve this 
problem. Although there are elements of the behavior 
of backflow which we do not presently understand, we 
come to the same conclusion as with X-shaped sources: 
the backflow model can account for other bridged distor- 
tions only in tandem with another hydrodynamic mech- 
anism. 

6. SUMMARY AND CONCLUSIONS 

We have conducted a series of three-dimensional hy- 
drodynamic simulations of light, hypersonic jets to study 
the viability of the backflow model for the formation of 
wings in X-shaped radio galaxies. The XRGs seem to be 
a population unto themselves, with characteristic envi- 
ronmental geometry, radio power, black hole mass, etc., 
and any successful model must account for these pecu- 
liarities. Our main results follow. 

1. The jets, back-flows, and lobes in our simulations 
are similar to those in the recent literature, giving 
us confldence in the usefulness of our simulations 
as probes of the backflow model. 

2. Wings in our models form in two stages: the es- 
tablishment of channels or proto-wings and then 
buoyant (usually subsonic) expansion. Speciflcally, 
our mod els corroborate the o verpressured cocoon 
model of iCapetti et al.l (|2002l ) early on, but as co- 
coons quickly come to pressure equilibrium with 
their surroundings, most of the wing length repre- 
sents subsonic expansion (Figure [3]). 

3. We have produced prominent wings by geometry 
and radio power alone, proving that the backflow 
model can, in principle, make X-shaped sources 
(Figure 12.21 and [TT]) . Both the atmosphere and jet 
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kinetic luminosity as a function of time are crucial 
to forming X-shaped sources. 

4. Long wings are produced in a relatively small por- 
tion of parameter space, requiring galaxies with 
high ellipticity, decaying jets, proper jet orienta- 
tion, and appropriate atmosphere size (Figures [7] 
and HI). 

5. The main challenges to the backflow model are the 
requirement for high ellipticity and subsonic wing 
growth. Adding additional physics is not obviously 
helpful. The backflow model seems to require an 
additional mechanism to make proto-wings which 
the backflow reinforces; in our models, we make 
these channels by the initial expansion of an over- 
pressured cocoon in an anisotropic environment, 
but this cannot explain every XRG. We cannot 
form new channels solely by deflecting back-flows. 

6. If the backflow model can overcome the issues 
noted above, it is a very strong candidate for ex- 
plaining X-shaped and other disturbed radio galax- 
ies. Our models naturally reproduce many of the 
characteristics of the XRG population (Table 2]) . 

There are several natural extensions of this work which 
promise to be fruitful. First, adding magnetic fields 
and investigating other potential coUimating mechanisms 
may rule out or boost the backflow model depending on 
the wing expansion speeds attained. Second, our mod- 
els rely on the formation of channels misaligned with the 
jets; these have been proposed in several other contexts 
as well to explain radio galaxy morphology. The origin of 
these channels is not known, and identifying and testing 
candidates would be important for wing formation mod- 
els. Third, more realistic jets could determine whether 
the backflow model can support both long wings and ac- 
tive FR II lobes with hot spots; the bending of the jets in 
particular is an important issue. Finally, models which 
produce the wings hydrodynamically but not explicitly 
by the deflection of bac kflow (e.g. the recently pro- 
posed stellar shell model iGop al-Krishna fc Wiital [2010( ) 
are worth exploring by making the atmospheres more 
complex in tandem with more realistic jets. 

The authors thank the anonymous referee for helpful 
and clarifying comments. We acknowledge support from 
Chandra grant GO011138A. 
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